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ABSTRACT 


This four- volume final report is concerned with Phase II of 
the DOE/JPL project s5 RF Model of the Distribution System As a 
Communication Channel." An earlier Phase I effort was concerned with 
the design, implementation, and verification of a computerized model 
for predicting the steady-state sinusoidal response of radial (tree) 
configured distribution feeders. That work demonstrated the feasi- 
bility and validity based on verification measurements made on a 
limited size portion of an actual live feeder. The Phase II effort 
is concerned with 1) extending the verification based on a greater 
variety of situations and network size, 2) extending the model capa- 
bilities for reverse directly* propagation, 3) investigating para- 
meter sensitivities, 4) Improving transformer models, and 5) investi- 
gating procedures/fixes for ameliorating propagation "trouble spots." 


PREFACE 


This volume contains more detailed (than Volume 1) accounts 
of the major tasks comprising this PHASE II effort. The reports 
comprising Task (1) and some of Task (2) (B) include verbatum 
copies of individual internal reports (by R. Wooding). 
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FOREWORD 


The detailed reports to follow for Subtasks ( 1 ) { A ) , (1)(B), and 
(1 ) (C) are copies of informal internal memos which were written at the 
completion of each task. Because of several contingencies which developed, 
these tasks were not written in chronological order of their numbering. 

The two principal contingencies requiring refinement of the model were: 

a) Discovery that original cable parameters were inaccurate 
(requiring development of a cable parameter program). 

b) Discovery of errors in line length specification using 
NMPC documentation. 

The first caused considerable delay and was the motivation for the tests, 
reported in Task (1)(C), on an electrically isolated 2-mile underground 
getaway cable. Similarly, discovering that line lengths as reported in 
utility documentation is not always accurate, especially for UG cables 
having much slower speeds of propagation, caused considerable delay. In 
fact, these two items introduced considerable confounding in early attempts 
to "unravel” the cause of observed discrepancies between predictions and 
measurements . 

For reasons mentioned in the discussions of subtasks ( 1 ) ( D) and 
(1)(E), planned work on these tasks was not accomplished. The discus- 
sions presented are verbatum from Volume 1. 

Finally, it is noted that during the process of making verifica- 
tion computations many software modifications were made to the DIFNAP 
system. Some of the most significant were: 

a) Adding a single phase ratio transformer capability 
(software & file) required for subtask ( 1 ) ( A) . 
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b) Improving the interactive communication logic 

c) Significant improvement to internal logic {for example, 
cumulative voltage transfer matrices) 

d) Simplification of distribution transformer logic and 

files for IDPL = 3. \ 

e) Significant reduction in record size for the "NT 

random binary file 

f) Similar reduction in the "DPU " record size. 

g) Including knowledge of actual frequency (in addition to code), 

passed down automatically through the "DPU „ 

This enables use of lumped parameter transformer models. 

h) Incorporating software for generating ZETA plotting. 

i) Incorporating a procedure for quickly introducing detailed 
outputing. 
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MODEL VERIFICATION - 
DENSELY LOADED FEEDER SECTIONS 


INTRODUCTION 

Distribution line carrier voltage propagation measurements were 
performed on a densely loaded overhead wire distribution circuit. The 
feeder section of interest is part of the 34556 feeder which originates 
from Niagara Mohawk Power Company's Grooms Rd. Substation in Clifton 
Park, N.Y. The nominal operating voltage of this feeder is 13.2 kv 
line to line. Figure 1 shows the general location of the densely loaded 
section in relation to the total feeder. Figure 2 is part of an opera- 
ting map provided by the utility showing the densely loaded section in 
more detail. This feeder section was selected for study because of the 
relatively large number of distribution transformers per unit length of 
ci rcui t. 

BACKGROUND 

A carrier voltage was placed on the feeder phase conductors at a 
capacitor bank location, shown at Point 1, Figure 2. Transmission was 
accomplished by connecting the transmitter as shown in Figure 3, thereby 
using the power factor correction capacitors to couple the carrier 
signals onto the phase conductors. The inductor was placed in the 
capacitor bank ground lead in order to prevent the bank from providing 
a large current sink to the carrier signals. The transmitter voltage 
was monitored at this point as shown in Figure 3. Transmission from the 
capacitor bank location was done as a matter of convenience. Except under 
unusual conditions, the inductor carrier voltage can be assumed to be 
very close to the voltage which appeared on the phase conductors. 
Transmission from the substation would have involved mounting separate 
capacitors and traps for each phase to determine the individual phase 
voltages impressed on the feeder section of interest. 
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FIGURE 5 


A coupling network with appropriate carrier frequency measurement 
equipment was connected near the end of the densely loaded section. 

This measurement apparatus, shown In Figure 4, was located at Point 2, 
Figure 2. 

A power factor correction capacitor bank was located at another 
terminal point of the densely loaded section. This bank was located 
at Point 3, Figure 2. An inductor was placed in the ground lead of 
this bank also for the reason mentioned above. The carrier frequency 
voltage across this Inductor was recorded during the measurements. 

The end result of the field measurements were voltage transfer 
ratios. Individual phase voltages at the transmission point were cal- 
culated by use of the measured transmit voltage at Point 1, Assuming 
the admittance of the power factor correction capacitors was much larger 
than the admittance seen looking out onto the feeder network at the 
capacitor bank location, the phase voltages at the capacitor bank are 
nearly equal to the Inductor voltage at the transmission point shown 
In Figure 3. Knowledge of these phase voltages at the transmission 
point provides a means to calculate individual voltage transfer ratios 
for each phase between Point 1 and Point 2 on Figure 2. Also, since the 
inductor voltage (a vector sum representative of the phase voltages at 
that point) was measured at the capacitor bank at Point 3, voltage 
transfer ratios between Point 1 and Point 3 were calculated from the 


measurements. 
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Figures 5, 6, and 7 show the comparisons of the measured and 
computed voltage transfer ratios between Points 1 and 2 for Phases A, 

B, and C, respectively. Figure 8 shows a comparison of measured and 
computed voltage transfer ratios between Point 1 and Point 3, the 
voltages being compared in Figure 8 are the Inductor voltages at the 
capacitor bank locations. (Point 1 and Point 3, Figure 2), 

CONCLUSION 

The model simulations show agreement of sufficient accuracy with 
the measured data for communication system engineering purposes. The 
model results are influenced somewhat by the loading level of the dis- 
tribution transformers. These simulations assumed a 20% load factor on 
the distribution transformers. This loading level was assumed after 
conferring with utility personnel, taking into consideration weather 
conditions and the time of day at which the measurements were made, 
with regard to system electrical loading. The power factor of the load 
on the distribution transformer secondaries was assumed to be .8 lagging. 
It should also be kept in mind the model assumed a time invariant network 
with respect to loading, while during the field test there probably was 
some variation in the feeder loading. 
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TASK (1 )(B) 


DIFNAP MODEL VERIFICATION 
FULL SIZE DISTRIBUTION FEEDER 


Introduction 

Distribution line carrier voltage propagation measurements 
were performed on feeder 34556 of the Niagara Mohawk Grooms Road 
Substation in Clifton Park, N.Y. The feeder operating voltage is 
13.2 kv line to line with single phase sections where the operating 
voltage is reduced to 4.8 kv line to line. These sections are shown 
in Figure 1. Simulations of these measurements were performed using 
the DIFNAP system programs in order to verify their utility in pre- 
dicting distribution line carrier propagation as an aid to communica- 
tion system engineering. 

Background 

Feeder 34556 is a radial distribution circuit. The total length 
of the main trunk circuit is about 10 miles. The following common 
distribution system features are found on 34556. 

1) Power factor correction capacitors 

2) Direct buried underground cables 

3) Ratio bank transformer's (13.2 kv - 4.8 kv) 

4) Long single phase branches 

5) "Densely" loaded residential areas 

6) 30 commercial & light industrial loads 

7) Several different overhead conductor geometries and 

conductor sizes. 
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The points labeled 1,2,3 & 4 on Figure 1 denote the location of 
30 capacitor banks. An inductor was placed in the ground lead of the 
capacitor bank's as shown in Figure 2. This was done in order to 
prevent the capacitors from providing a large current sink to the 
carrier signals. During the measurements, the carrier voltage across 
the inductor was monitored and recorded to give an indication of the 
vector sum of the phase voltages. Sufficient precautions were taken 
to prevent 60 Hz. voltages from affecting the carrier voltage measure- 
ments . 

Individual phase voltages were monitored at point (2), Figure 1. 
At this point, access was gained to the individual overhead phase 
conductors. Within a few hundred feet of this pole is a cable riser 
providing service for the Country Knolls West residential area which 
is served by direct buried underground cable. Extensive measurements 
and DIFNAP simulations on this cable network were reviewed in a 
previous report. The coupling network employed at Carlton Rd. (point 
2, Figure 1) is shown in Figure 3. 

During the field measurement, three distinct conditions existed 
on the feeder which were treated by simulation. These conditions were: 

Condi tion 1 . No coupling network connected at Carlton Rd. 
Inductor voltages were monitored at Points (?.) and (3) shown 
on Figure 1. Transmit voltages were monitored on Phase A 
at the substation as shown in Figure 4. Phase voltages were 
assumed equal at the substation for modeling purposes. 

Condition 2 . Coupling network connected at Carlton Rd., but 
with tuning inductors disconnected. Individual phase voltages 
were monitored at Carlton Rd. Inductor and transmit voltages 
monitored as in Condition 1, 
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Condition 3 . Coupling network" connected with tuning inductors 

in place at Carlton Rd. The same quantities mentioned in 

Condition 2 were again recorded. 

Simulations were performed at the following frequencies: 

5010 Hz 
6990 Hz 
8130 Hz 
9510 Hz 

Corresponding to measurements performed on the feeder. Table 1 
provides the diagonal elements of the admittance matrices representing 
the coupling/tuning network loading under Conditions 2 and 3. The 
off-diagonal elements of these matrices were assumed to be negligible. 

Resul ts 

The first situation to be considered will be the Condition 1 out 
lined above. Under this condition, the network is essentially "normal 
Figures 5 and 6 show a comparison of the simulation results and field 
measurements of the inductor >jltages at points 1 and 3, respectively. 
The inductor voltage was arrived at by making use of the admittance 
matrix (Y) of the inductor-capacitor network. The current vector 
I is, I = Y Ep, where is the phase voltage vector provided by 

the program. The sum of the currents in the vector I was then used to 
calculate the voltage across the inductor. The voltage transfer ratio 
to the capacitor bank at Point 1 is much higher and more uniform 
across the frequency range than that at Point 3. 

Under Condition 2, the coupling network was now attached to the 
phase conductors at Point 2 on Carlton Rd. This coupling network. 
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diagrammed in Figure 3, represented a significant loading effect on 
the network. Table 1 shows that without the tuning inductors attached, 
the coupling network presents an inductive load. Figures 7 & 8 show 
comparisons ojf measured & computed'results" at Points (V) & (3). Wfth 
the coupling network now attached, it was possible to measure the 
individual phase voltages. Figures 9, 10 and 11 show comparisons 
between measured and computed voltages on Phases A, B, and C of 
Carlton Road under Condition 2. 

Under Condition 3, the tuning Inductors were attached to the 
coupling/network. These inductors were then varied independently in 
order to maximize the voltage transfer from the substation to the 
capacitor bank at Point 3, Figure 1. Table 1 shows that Phases A and 
C were loaded capacitively, while Phase B was loaded inductively at 
Carlton Rd. at 5010 Hz. The tuning elements then remained fixed as 
frequency was increased. Figures 12 and 13 show the results obtained 
for the inductor voltages at Points 1 and 3, respectively . Figures 14, 
15 and 16 show the phase voltage results under Condition 3 on Phases 
A, B, and C, respectively. 

Conclusions 

The simulations performed on the 34556 distribution feeder and 
subsequent comparisons to measured results were subject to several 
assumptions. These being: 

1) The network configuration was specified accurately by 
utility personnel and on system maps. 

2) Conductor lengths specified on utility map were accurate 
for engineering purposes. 
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3) The assumption that each transformer was loaded to 20% 

of full kva capacity with a .8 power factor lagging load 
was a valid representation for distribution transformer 
loading. 

Under this assumption, the total connected load kva 
employed in the model was 3200 kva. If the load factor 
mentioned above was assumed. The measurements were per- 
formed during off-peak loading hours for the most part under 
moderate climatic conditions with respect to heating and 
cooling loads. 

4) The coupling and tuning network components were known with 
sufficient accuracy for model purposes, and loss components 
present In these circuit elements were negligible at the 
frequency of interest. 

5) Errors In field voltage measurements and data reduction 
were minor. 

Since the voltage transfer ratio of the feeder is sensitive to the 
loading level, it would be advisable for engineering purposes to study 
the feeder response parametrically as a function of loading level in 
order to bracket system performance limits. 
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TASK (1)(C) 

Underground-Overhead Transi ti on 


Introduction 

As previously mentioned in the Foreword, two significant con- 
tingencies were encountered during the first attempts at verification. 
Both of these occurred with "subnetworks" involving underground cables. 
After considerable delay and extra experimentation, these contingencies 
were found to be caused by; 

a) Inaccurate cable parameters 

b) Errors in determining actual line lengths from utility 
documentation. 

The first arose in attempting to adapt the Bonneville EMTP 
program to compute parameters for a cable configuration which was not 
explicitly available in its repertoire. To overcome this difficulty, 

a cable parameters proqram was developed to handle single con- 
ductor concentric, non ferrite sheathed cable. Also, it was discovered 
by making tests on samples, that the dielectric had a much higher 
dielectric constant than was originally assumed. Finally, it was 
discovered that line lengths given on utility documentation often did 
not take into account extra lengths for short laterals to vaults and 
transformers, and that often significant cable length was left in terms 
of a coiled end (for future splicing). 

The contingencies resulted in some laboratory tests on cable 
samples, development of the computer program, and a special field test 
on a 2-mile underground cable which was electrically isolated and made 
available by NMPG for testing. 


The results of three separate field tests related to underground 


and underground-overhead transitions are included herein. 


Model Verification G.R.S.S. U-G Getaway 


INTRODUCTION 

This report describes the procedures undertaken to verify a 
computer model of measurements taken on an underground trunk line 
feeder cable. The 30 cable system treated in this verification effort 
serves one of the several distribution feeders in place at the Grooms 
Road Substation of Niagara Mohawk Power Co. 

BACKGROUND 

A diagram of this feeder configuration is shown in Fig. 1. The 
trunk line underground feeder leaves G.R.S.S. (Grooms Road Substation) 
and continues underground for a distance of 1.98 miles, after which 
the cable surfaces at a riser pole to feed a hybrid overhead-underground 
distribution system. The trunkline feeder cable was isolated at each 
end for the purposes of measurement and verification. When isolated, 
the trunk line cable therefore constitutes an unbranched 30 trans- 
mission line with various possible terminations. 

The main conductor size of the cable is 750,000 CM and the 
material is aluminum. The neutral conductor is composed of 21 
#14 copper wires wrapped concentrically over the main insulation 
jacket as shown in Figure 2. The main conductor is covered by a semi- 
conducting jacket and the main insulation jacket is covered by a semi- 
conducting jacket as shown in Figure 2. These jackets are present 


1-38 


ORIGINAL PAGi* 

OF POOR QUALITY 



SUBSTATION 


FIGURE 1 


POLYETHYLENE JACKET 


NEUTRAL WIRES 


EMI- CONDUCTING JACKET 



MAIN CONDUCTOR 


MAIN INSULATION 


FIGURE 2 
1-39 



for electric field grading purposes* The neutral wires are covered by 
a polyethylene jacket. The main insulation system is composed of 
"Kerite", whose relative permi ttivity is estimated to be about 3,0 

These conductors are buried to constitute a 30 circuit in equi- 
lateral fashion as shown in Figure 3. The burial depth is a minimum 
of 42 Inches. A solid soft drawn #2 copper wire is buried along with 
the cables to serve as a safety ground. The cable neutral wires are 
bonded to the copper wires at any splice point, which are a maximum of 
1000’ apart. 

In constructing the computer model of this cable, several factors 
which affected the model results were not known with a high degree of 
certainty. Among these factors were: 

1) Earth resistivity 

2) Main conductor and neutral wire temperature 

3) Physical spacing of conductors 

4) Exact dielectric constant of insulation system 

5) Actual cable length 

Efforts were made to employ the cable parameters portion of the 
E.M.T. program of the Bonneville Power Administration in order to obtain 
the necessary transmission line parameter data base files. The parti- 
cular configuration of neutral wires found in these concentric neutrcl 
cables was not directly amenable for calculation by use of the EMTP. 

A cable parameter program was developed to handle modeling of 
this type of cable construction. 
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To verify the underground cable model, calculations were made 
for the 30 driving point admittance matrix corresponding to that 
measured at GRSS with the various terminations applied to the end 
of the Isolated underground "getaway" trunk cable shown in Figure 1. 

The cable parameters and the driving point admittance matrix were 
calculated for the following frequencies: 

6440 Hz 
8130 Hz 
9510 Hz 
12010 Hz 
17010 Hz 
20010 Hz 
21980 Hz 

The simulations made at these frequencies were compared with results 
from the measurements at similar frequencies. According to the net- 
work maps provided by N.M.P., the cable length was taken to be 3189M. 

In general, the measurement of the main diagonal elements of the admittance 
matrix for a general 3 terminal network was undertaken as shown in 
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Figure 4, where * y— - n 

The cable network model was first configured to simulate the 
situation where the terminal end of the cable was open-circuited. 

The results of this computation as a function of frequency are shown 
in Figure 5 along with the results of the field measurements. Note 
that the computed admittance magnitude null occurs below that of the 
measured null. A 5 percent error was assumed in the estimation of the 
length of the network (5% too long), and the admittance values recom- 
puted with the result shown in Figure 6. In Figure 6, it can be seen 
that this length modification brings the computed and measured results 
into good agreement. Due to the uncertainty in the precise length of 
the cable, and the other important factors mentioned earlier, this re- 
computation, for purposes of the model verification, is not unreasonable. 

Measurements were made to determine the value of the main diagonal 
of the admittance matrix as a function of frequency with a resistive 
termination network placed at the end of the cables as shown in Figure 8. 
The computed and measured response for this situation is shov/n in 
Figure 7. By comparison with Figure 5, it can be seen that this termina- 
tion eliminated some standing wave effects seen in the unterminated case. 

Computations performed by the FEEDPUS3 program show the charac- 
teristic, impedance of the cable under study to be about 15.5Q. A 15n 
termination was selected. 

In the general case, the measurement of the off-diagonal term 
of the admittance matrix is shown in Figure 9. 
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The computations performed by the NTWKERH1 program show the 
magnitude of the off-diagonal terms to be on the order of -40 db,, 
smaller than those of the main diagonal. The measurements, in general, 
show this term to be only about -20 db smaller than the main diagonal 
term at first glance. Figure 10 shows the measurements of an off- 
diagonal term for the terminated cable situation. Notice that the 
phase measure shows a continuous accumulation of the phase difference 
between the referenced voltage and the measured phase current. It should 
also be noted that the two 360° phase rolls occur, on the average, at 
about 17.2 kHz intervals. Calculations from FEEDPUS3 program show 
the velocity of propagation to be about .35 the speed of free space 
propagation for this cable. 

From this data one may calculate the 2-way travel distance 
V7 ~ 2 ~ kH z = 5 - 8139 sec x 10 “? D=VT=.35( 3 x 10 8 ) (5.81 34 x 10" 5 sec) = 6104M. 

Now if the current measured in finding the off-diagonal term is the 
return currenttraveling out thru one phase conductor and returning thru the 
other, as shown in Figure 11, rather than the through ground path, the figure 
obtained above, when divided by 2, should equal the one way distance of the 
cable. 

So, 

= 3052M 
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FIGURE 10 






When compared to the assumed 301 1 M . cable length used in the main 
diagonal measurements, this compares quite favorably. Therefore, if 
this return current is present, the difference between calculated and 
measured off-diagonal terms is not surprising. In order to investigate 
the general possibility of simulating the transport delay situation 
described above, a delta load was provided as a termination using the 
computer model . The results obtained from frequency response calcula- 
tions using the model showed that this general response pattern could be 
simulated. No attempt was made to simulate the response exactly due to 
lack of information regarding the exact coupling mechanism. Also, since 
the computed value of the off-diagonal terms are small relative to the 
main diagonal terms, the transport mechanism noted above, whatever that 
may be, could easily mask the off-diagonal measurement. This inability 
to measure and verify the computation of the off-diagonal term of the 
admittance matrix is not serious since the degree of coupling between the 
cables is so low initially that the off-diagonal terms will probably make 
little difference In general network calculations. 

Measurements of the off-diagonal term were taken also for the 
situation where the cable was open-circuited at the endpoint. The results 
of these measurements are shown in Fig. 12. By examination of the phase 
angle plot of the off-diagonal term, it can be seen that the transport 
mechanism has not completely masked out the mutual coupling between the 
phase cables. 

SUMMARY 

In summary, the difficulty in verification of the off-diagonal 
terms of the cable admittance matrix is not seen to be a detriment to 
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the overall modeling effort since the off-diagonal terms are relatively 
small with respect to the main diagonal of the admittance matrix. Also, 
it would not seem to be likely that this transport mechanism would 
exhibit itself under actual network condition, since any large loads 
{P.F.C. Capacitors, etc.) which provide coupling between phase conductors, 
would tend to dominate the off-diagonal component of the admittance matrix 
due to Intercable coupling. 


MODEL VERIFICATION OF MEASUREMENTS PERFORMED 
ON COUNTRY KNOLLS SOUTH U.G. NETWORK 


INTRODUCTION 


BACKGROUND 


This report describes procedures undertaken to verify 
field measurements made to determine the driving point 
admittance of an underground tree network at several 
discrete frequencies using the DIFNAP computer 
program system. 

Field measurements were performed in order to measure 
the elements of tins admittance matrix, at various fre- 
quencies, of an entire underground tree type distribu- 
tion network. 

In general, the network consists of a large trunk 
cable with a number of single phase branches implemented 
with smaller conductor cable. The main trunk cable is 
about 1600 meters long and is composed of 500 KCMIL con- 
centric neutral cable. The single, phase branch circuits 
are composed of #2 A.W.G. concentric neutral cable. 

In order to measure the elements of the main diagonal 
of the admittance matrix, the general configuration shown 
in Figure 1 was employed. A voltage was placed on one of 
the phases while the other two phases should be at refer- 
ence potential. The phasor ratio of the voltage and current 
on the excited phase yields the value of the main diagonal 
element. In practice, the voltages on the non excited 
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phases were not able to be lowered sufficiently close to 
reference potential. This difficulty arose due to the fact 
that the distribution feeder circuits were measured while 
energized, eliminating the use of direct grounding. During 
the initial field measurements, a series L-C tuned circuit 
was employed to provide a low impedance path to ground on 
the unexcited phases for the carrier signal. This technique 
provided a maximum 20 db differential between the voltage on 
the excited phase and that of the "grounded" phases, with the 
norm being closer to 15 db differential. In order to provide 
a more suitable "ground" on the unexcited phases, a voltage 
cancellation technique was used on the unexcited phases. 

During follow-up field measurements, this technique forced 
the voltage differential Into the neighborhood of 45 db. 

The measurement of the off-diagonal terms was undertaken 
as shown in Figure 2. 

Figure 3 shows the results of the field measurements of the 
main diagonal element of the admittance matrix for Phase A 
of the underground network in addition to the values computed 
by use of the NTWKERH1 computer program. Values are shown 
for measurement and computations between the frequencies 
of 5 kHz and 25 kHz. 

Field test results obtained using both nulling techni- 
ques are employed. The results from the first 
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set of field measurements are used even though they are less 
desirable than the results of the second set, since the 
second series of measurements measured the following fre- 
quencies between 5 and 25 kHz: 


5010 

HZ. 

6990 

HZ. 

9510 

HZ: 

21010 

HZ. 

23010 

HZ. 

25010 

HZ. 


These frequencies were not sufficient to fully reflect 
resonance effects in the cable network so results from the 
first series of tests were included. Frequencies measured 
in common on both tests show good agreement, taking into 
account the dynamic nature of the network with respect to 
relative loading levels. 

Figures 4 and 5 show the results of the measurements 
and computer simulations of the main diagonal elements of 
the admittance matrix for Phases B & C of the underground 
cable network. 

From inspection of Figures ,3, 4, and 5, it can be seen 
that there is good agreement between the frequencies at 
which the relative minima and maxima of magnitude and 
phase occur. 
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A distribution transformer loading level of 4Q5» of full load 
rating was selected on the basis of the time of day at which the 
tests were performed and weather conditions. 

It can be seen from the results shown in Figures 3, 4, and 5 
that the model is quite useful in simulating the response of the 
U.G. cable network considering the uncertainties noted above and the 
accuracy needed for communication system design. 

Sufficient data was taken to compute the off diagonal terms 
of the admittance during the field tests. The results of the field 
test show the magnitude of the off diagonal terms to be about 30 db 
smaller than those on the main diagonal. Results of the computer 
model simulation show these off diagonal terms to be on the order 
of 60 db smaller than those on the main diagonal. These results 
reflect the high degree of shielding between the phase cables. 
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MODEL VERIFICATION OF FIELD MEASUREMENT 
PERFORMED ON CARLTON RD. O.H. - U.6. NETWORK 


INTRODUCTION 

Sufficient field data was gathered to provide a measure of the 
complete admittance matrix of the overhead-underground network origi- 
nating approximately at the junction of Carlton Rd. and the Jonesville- 
Elnora Rd. in Clifton Park, on the NMPC 34556 feeder. 

In approximate terms, the network consists of 650 meters of 30 
open wire primary feeding a large underground trunk cable of 500 meters 
length. Several 10 feeder cables are branched from this trunk and are 
about 600 meters in length. The trunk cable consists of 500 KCM 
concentric neutral cable. The feeder cables are #2 A.W.G. size, with 
concentric neutrals. 

A description of the definition of the elements of the admittance 
matrix is given in two earlier reports describing verification efforts 
on underground cable networks. 

During the admittance matrix measurements, each phase was excited 
individually, while the other two phases were theoretically brought to 
signal reference potential. While transmission occurred on each phase, 
secondary line to neutral voltage was measured on a distribution trans- 
former on a single phase branch. These secondary DT voltages were used 
to infer the primary voltage on the cable by adding a turns ratio correc 
tion to this voltage (36 db). Using these estimated values of primary 
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voltage, a measure of the voltage transfer ratio of O.H. - U.G. network 
could be obtained. The admittance matrix was then calculated using the 
NTWKERH1 program while the voltage transfer was calculated using the 
NTWKANH1 program. 

Field measurements were also performed to collect data on the above 
mentioned quantities when a reactive load {composed of the test van 
coupling network) was placed on the end of the trunk cable. The reactive 
load on each phase of the trunk cable was varied in order to maximize the 
voltage transfer outbound between the substation and a capacitor bank 
location (capacitor bank at point K, Figure 1). The reactive loads were 
varied as a crude experiment to see if voltage transfer between the sub- 
station and point K could be influenced by the reactive load at the end 
of the trunk cable. Simulations were performed on the network for this 
condition also. 

All D.T.'s on the network were assumed to be loaded to 40% capacity 
with a unity power factor load. This loading level was selected with due 
consideration of time of day and weather conditions as influences on 
electricity usage. 

Results 

Figures 2, 3, and 4 are comparisons of the measured and computed 

values of the main diagonal of the admittance matrix when the network was 

not modified. The frequencies at which measurements were taken were: 

5010 Hz 
6990 Hz 
9510 Hz 
21090 Hz 
23070 Hz 
25510 Hz 
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N.M.P.C. 34556 FEEDER 
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COMPUTED MAGNITUDE 



Simulations and measurements were performed for the off-diagonal 
elements as well. In this case, a 3x3 admittance matrix is involved. 

The Y-| 2 » Y^> Y 23 elements were examined. Figures 5, 6 , and 7 show 
model results and corresponding measurements . 

Figures 8 and 9 show results of voltage transfer measurements 
and modeling. As mentioned above, carrier signal voltages were moni- 
tored on the secondary of a distribution transformer. These voltages 
were measured line to neutral. In order to obtain a rough idea of the 
voltage on the primary conductor, a correction for the transformer turns 
ratio was added (36 db) to this secondary voltage. The NTWKANH1 program 
was then used to perform the voltage transfer simulation. The D.T. was 
located on a single phase branch being fed by Phase One. Figure 8 
compares measured and computed results of the voltage transfer ratio 
while transmission occurred on Phase 1. Figure 9 shows the same com- 
parison while transmitting on Phase 2. It should be kept: in mind that 
no attempt was made to account for the voltage transfer characteristics 
of the D.T. at the various frequencies of interest, only a turns ratio 
correction was used. 

Also, as mentioned earlier, a reactive load was placed on the end 
of the trunk cable. Simulations and measurements were performed for this 
condition at the following frequencies: 


5010 Hz 
6990 Hz 
9510 Hz 
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FIGURE 


Figures 10, 11 and 12 compare the measured and modeled main 
diagonal elements of the admittance matrix for this modified network 
condition. The load on each phase was different and was chosen in 
order to optimize phase voltages further out on the feeder by trial 
and error methods. 

Figures 13, 14 and 15 compare the results of measurement and 
simulation for the Y^. Y^, and ^22 elements, respectively, for this 
modified network condition. 

Again, measurements were taken on the D.T. secondary while the 
network was In a modified condition. Figure 16 shows the voltage trans- 
fer ratio comparison of the computer simulation with the inferred 
primary voltages from the D.T. while transmitting on Phase A. 

Figure 17 shows a similar comparison for transmission to Phase B. 


Conclusions : 

The "Hybrid" overhead-underground network lent itself well to 
simulation by the NTWKERH1 and NTWKANHl programs. As opposed to a pure 
cable network, the relatively short overhead section provided a substan- 
tial amount of ii ^erphase coupling. (Studies of U.G. cables indicate 
little mutual coupling). This is illustrated by the fact that reasonable 
signal levels could be obtained on the Phase 1 *,T. while transmitting 
on the other phases during the admittance matrix measurement. 
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The voltage transfer ratio comparison showed good agreement, 
considering the crude method used to infer the primary voltage at the 
D.T. during the measurements. 

Also, it can be seen that it is quite easy to simulate the effects 
of the reactive load modification of the network with good accuracy. 

Again, it should be noted that various uncertainties are present 
in the network model such as: 

1) Network length uncertainty 

2) Conductor temperature 

3) System loading (size and location) 

4) Earth conductivity 

Nevertheless, the results of this simulation would be adequate for 
communication system engineering. 
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TASK (1 )(D) 

Meter Terminals on Secondary* 


The purpose of this subtask was to verify the models ability 
to predict voltages at the "meter terminals" on the secondary service 
of distribution transformers. Inasmuch as it was not planned as part 
of this PHASE II effort to have a distribution transformer analytic 
model implemented, the equivalent prediction is presently performed by 
a "two-step" process. The first step is to predict propagation on the 
feeder with distribution transformers (DT) represented by a phasor 
admittance Toad connected to the proper phase and representing the DT 
primary driving point admittance associated with the specified DT and 
connected secondary load. Then having the predicted feeder phase 
voltages, an off-line calculation or response curve look up could be 
used to predict the DT secondary terminal voltage. The voltage propa- • 
gation down to the meter terminals could also be done off-line using 
electrically short (primarily a lumped series impedance) secondary 
circuit representation. 

Originally it was hoped to find a DT on NMPC feeder 34556 which 
had a single customer at the end of a long secondary service. To 
verify with a "multi drop" secondary would have required simultaneous 
measurement of each customer load-current and voltage which was not 
considered feasible. However, no long single customer service was 
found. Therefore, no uniquely applicable test and verification was 
conducted. 

However, several related activities and characteristics have some 
positive bearing on this objective. First, as part of TASK 1C, propa- 
gations were made from a location on the feeder (removed from the DT 


♦This discussion is verbatum from Volume 1 
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location) to the secondary terminals of various DT's, and the results 
were supportive. Secondly, as part of TASK 4 DT modeling extensive 
response measurements were made that show, for the lower frequency 
range (say less than 25 kHz) the primary to secondary transfer is 
primarily the turns ratio effect with & few DB accounting for reactance 
drop caused by secondary load. Finally, any voltage drops on the 
secondary circuit would lepend on an estimate for the simultaneous 
loads for each customer. Since the secondary circuits are electrically 
short, the computation, given the loads, presents no conceptual problem. 

It is believed, therefore, that compared to predicting propagation 
on the feeder, analysis of the DT primary to secondary (terminals) 
transfer and prediction on the secondary terminals present no conceptual 
difficulty. 
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TASK (1)(E) 

Extended Frequency Range* 


The objective of this subtask was to extend the frequency range 
of verification to satisfy the objective of PHASE I. Because of 
technical considerations mentioned in 1) and 2) below, it was deemed 
higher priority to verify first the models capability at the lower 
frequencies. Hence, the other tasks associated with verifying the 
model were given scheduling priority. For the reasons cited below, 
the work called for by this task was not accomplished. 

1) Uncertainties in Specifying Line Lengths 

As work progressed to the point of verifying the model at 
lower frequencies, it was realized that apriori estimates {based on 
utility records) of line lengths are not always reliable. This aspect 
required a non trivial effort to understand the reason for lack of 
verification in some of the trials. Because of wavelength considera- 
tions, errors in assigning (to the model) line lengths will become 
"damaging" with higher frequency. Based on our experience, this 
limit may be 50-100 kHz. 

2) Transformer Response 

As the transformer modelling- and measurements progressed, it 
appeared more and more promising that with a simple extension (to the 
60 Hz R-L-M two winding transformer model, a lumped parameter model 
could be used to predict transformer response up to 25-50 kHz and which 
was not too sensitive to variations in construction, i.e., with the simple 
extension the similarity of 60 Hz performance could be extended to 25- 
BO kHz. (Characteristics do vary with KVA). Beyond this frequency 
range, variations in response at PLC due to differences In manufacturing 
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construction may become important. If this is true, than either a 
much more detailed inventory of transformer characteristics and identi- 
fication of installed transformers is required, or additional modeling 
errors accepted. 

3) Resource and Time Limitations 

During the process of verifying the model, several contin- 
gencies were encountered. One involving line length uncertainties was 
mentioned above. The most significant contingency involved verifica- 
tion with circuits involving underground cables. Identifying the source 
of the problem, arranging for appropriate teste (see Volume II, Task 
(1)(C), writing a new cable parameter program involved considerable 
effort and delay. Another contingency involved making admittance 
measurements of the Grooms Road Substation. Special measurement equip- 
ment had to be designed and built and special facilities installed by 
Niagara Mohawk Power. 

As a result, resources (and time) were not available to proceed 
with higher frequencies. 

4) Vendor Trends to Lower Frequencies 

Supporting the above stated priority is what appeared to be 
a trend by virtually all parties to gravitate to the lower frequencies 
for distribution feeder communication. 
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Summary of Verification Activities 

As a result of a considerable experience gained in verifying the 
model, the following caveats are offered:- 

1) Line length uncertainties, especially for underground 
sections, can be significant such that, typically, for 
frequencies greater than 25-50 kHz, these lengths should 
be correctly defined. Even at lower frequencies, point 
estimates of a given frequency should be supplemented by 
a frequency response in the neighborhood to reveal if in 
the vicinity of a resonance which could be sensitive to 
line length errors. 

The above mentioned frequency limit is also applicable if 
the simple lumped parameter transformer model is to be 
used. 

2) Line loading can have a significant effect (in reality and 
model) so it is recommended that this effect be investiga- 
ted parametrically for "worst case situations". 

3) Inbound calculations which usually involve representation 
of a substation (admittance load) will be sensitive to 
estimation errors in this admittance representation unless 
the design (trapping, etc.) negates the effect of these 
errors . 

Except for item 1) the above represent reality and as such are not a 
problem in the modeling. Hence, it is recommended that appropriate 
parameter variations (in relative and specific loading) be conducted 
to determine worst case effects. Item 1) suggests that in addition to 
using utility documentation, appropriate questions be asked of the 
utility operations to uncover possible practices which would in fact 
yield line lengths different (usually longer) than documented. 
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Having these caveats in mind, it is our opinion that the model 
predicts quite well and should be useful for communication design, 
and that if a continued need exists for this kind of tool that it would 
be justified to further enhance its capability and convenience (see 
section containing recommendations in Volume 1). 


1-38 



TASK 2 


EXTENSION, IMPLEMENTATION & VERIFICATION OF 
PERTURBATION THEORY 


2-1 


TASK 2 


TABLE OF CONTENTS 


Page 


FOREWARD 

2-1 

TASK 2(-) 

2-2 

TASK 2(A) 

2-8 


TASK 2(B) 


2-9 


CONTRIBUTORS 


TASK 2 (- } 
TASK 2(A) 
TASK 2(B) 


J. T. Gajjar 

R. Wooding, J. T. Gajjar 

R. Wooding, R. W. Rankin, M. J. Luka 
G. W. Kinloch, W. C. Hughes 


2-iii 



FOREWORD 


l 


2-1 


The work statement specified three types of activities 
associated with this task, i.e,, 

TASK 2(-) Extension and Implementation of 
Perturbation Theory 

TASK 2(A) Use Perturbation to Calculate Inbound 
Path Propagation 

Task 2(B) Measure Inbound Path Propagation 

The work associated with the first activity is reported herein 
supplemented by the material contained in Appendix 1 and 2 in 
Volume 3, The work and results associated with Task 2(A) and 

2(B) are contained herein. 
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TASK 2(-) EXTENSION AND IMPLEMENTATION OF PERTURBATION THEORY 


The DIFNAP system of programs is a useful tool for the analysis 
of signal flow over power distribution systems configured in the form 
of distributed, multi -conductor trees. The analysis procedure used 
in these programs is a tv.*!'}- step procedure. In the first step, starting 
with the extremities of the network, input admittances are computed" 
for each branch of the tree. Voltage transfer ratios across each 
section of the tree are also computed and saved. The process culminates 
in the computation of the input admittance to the entire network. In 
the second step, the input voltages are determined using the input admit- 
tance to the network computed in the first step, and the source para- 
meters. Voltages at various points of the network can be calculated 
using the input voltages and the voltage transfer matrice computed and 
saved in the first step. 

By computing and saving the admittance parameters of each section 
during the first step, currents and signal power at various points in 
the network can also be computed. 

It is often of interest to determine the effects of changes in 
load values on the voltages in the network. The use of the DIFNAP 
program would normally entail recomputing the entire network with both 
steps. 

In order to permit more efficient computations of effects of load 
variation, during the PHASE I model development, a technique 
called the "perturbation theory" was developed. This technique uses 
the application of the reciprocity theorem, as extended to multi- 
conductor systems, to the adjoint of the tree structured mul ti conductor 
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distribution feeder network. 

As a consequence, voltage variation at any point on the network, 
resulting from the variation of a single load, can be computed by 
exercising the second step of the procedure only. It was determined 
that the inbound computation of network input admittances and re- 
computation of voltage ratios is not necessary. Although we have 
chosen to call the process application of perturbation theory, since 
we are dealing with changes in voltage resulting from a change in a 
load value, the computation is in no way an approximation. In fact, 
the procedure is computationally exact, even for large variations in 
load. It can, for instance, be used for computing the effects of 
faults, both open and short, on signal propagation in the model. 

Further details of the mathematical development of the theory and 

(121 

algorithm can be found in the Sponsor Reports' * ' submitted in con- 
junction with the Phase I report. A tutorial on this subject, written 
by R. C. Rustay, has been included as Appendix 1 in Volume 3. 

The theory and algorithms were developed during the first phase 
of the model development. Time and resources did not permit us to 
incorporate the algorithms into the model during that phase. 

During the INTERPHASE period following the completion of PHASE I 
and before the beginning of PHASE II of the project, we realized the 
importance of an ability to make inbound signal propagation computations 
from points on the network to a receiver. It was also deemed essential 
that the location of the receiver be not necessarily the same as that 
of the source in the outbound computation. 
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One possible approach to this problem was to create a different 
network with each load point of the original network being the source 
point of the new network. By using both of the steps of the DIFNAP 
program procedures outlined above, the circuit can be analyzed and the 
inbound signal at the receiver located at any point on the network 
could be estimated. For a large number of load points, however, this 
procedure becomes both cumbersome and computationally uneconomical. 

The other approach is to make direct use of the reciprocal proper- 
ties of bilinear, passive networks. In this approach one uses transpose 
matrix relationships and computations made in the outbound path to 
determine the reverse path computations. The application of a source 
at a load point for the purpose of transmitting information, usually 
also results in a change in admittance at the load point. Thus, in 
addition to the application of reciprocity, a technique such as the 
perturbation procedure outlined above, has to be used to make the 
computations. This approach, however, suffers from the limitations 
that we can only make computations for the receiver located at the 
outbound source of the network. 

From the mathematical theory developed, for the perturbation 
approach during Phase I, we recognized the potential for extending 
this approach to the solution of the reverse path problem. This 
extension was undertaken and completed during the INTERPHA5E period 
and the perturbation algorithms together with those for the reverse 
path were implemented into the computer model. Just prior to the 
beginning of the Phase II effort, we also developed a program which 
permits the network reversal from any point on the network and using 
it we were able to prove out the procedure implemented for perturbation 


and reverse path propagation. 

The perturbation and reverse path procedures are implemented so 
that they can be entered after the initial processing of the second 
step for at least one outbound computation. At the end of each forward 
processing, the interactive program asks the user if a perturbation 
on reverse path computation is desired. 

When the reverse path/perturbation option is chosen, the user 
next enters the node number of the node at which the load variation 
takes place or the reverse path source is applied. After verifying 
this, the program internally checks to determine if the outbound path 
to this point has been computed. If not, the program proceeds to 
make the outbound computation to that point. 

The user is then asked to choose between a reverse path or per- 
turbation computation. If the reverse path computation is chosen, the 
program asks him to define the reverse path source as a Norton's 
equivalent source. If the perturbation computation is chosen, the 
user defines the change in load. Following the verification of the 
source, or load change, the user is asked to define the end of the 
section where the result is desired. The program then computes and 
displays the outbound voltages, change in the outbound voltage due 
to the load perturbation or the reverse path source, and the algebraic 
sum of the two. 

The user can then choose another point at which to evaluate the 
effect of the network change. If not, he can choose another point 
at which to insert a reverse path source or load perturbation. If 
he chooses neither, the program exits from the perturbation mode and 
asks the user if he would like to compute the outbound network with 
another source. 
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The user interaction of the analysis program (step 2) with the 
reverse path option is demonstrated in the printout shown in 
Figure 2^1_. The perturbation interaction is similar except that no 
current source input is involved. 

Operationally, the perturbation/reverse path section does 
not alter any of the information computed during the outbound path 
computation, as it does not involve the output macro in its execution. 
Thus, in repeated applications of the perturbation procedure, consistent 
results are obtained after several different applications. 

The perturbation/reverse path approach permits computation of 
the effects of introducing a source and/or a load admittance change 
at any point in the network. The point at which the result is deter- 
mined can be in another branch of the network or can be in an out- 
bound direction from the point at which the source/load variation 
are applied. Thus, point-to-point signal propagation and the 
solution of partial networks can also be accomplished by this 
approach. 
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Figure 2-1 
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Figure 2-1 
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TASK 2(A) PERTURBATION CALCULATION OF INBOUND PROPAGATION 


This task was concerned with using the implemented perturbation 
technique to predict inbound/reverss propagation. What actually was 
done was to predict inbound propagation using both perturbation techni- 
ques and an alternate procedure involving a "network reversal" operation* 
followed by a normal outbound prediction from the "new" source to, say, 
the substation. The purpose of using this alternate procedure was first 
to give an "exact" numerical check on the perturbation arithmetic imple- 
mentation and also as a more convenient procedure for the general veri- 
fication of the model. The result of this model verification is contained 
in the following Task (2) (B) Section. 

Since the arithmetic results using the perturbation and the 
reversed network procedures gave exactly (except for negligible 
differences due to numerical round off considerations) the same results 
in all cases, specific numerical results are not presented here. Also, 
the same procedures were used on other synthetic networks designed to 
test various topological considerations (conceptually three different 
topological classes exist, each of which requires somewhat different 
perturbation logic. Also within each of these classes, various other 
situations, such as transitions, were tested). As a result of these 
tests, it is believed that the implemented perturbation logic is 
validated. 


*A main program "NETREVSI" was implemented during the INTERPHASE interval 
and which can automatically generate from a given network file a new 
network file having a specified remote node as a new source. 
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MODEL VERIFICATION - INBOUND PATH PROPAGATION 


INTRODUCTION 

Inbound signal propagation measurements were performed on a 
13.2 kv distribution feeder. This feeder being the 34556 feeder on 
Niagara Mohawk's Grooms Road Substation in Clifton Park, N.Y. Inbound 
transmissions were made at a point about 7.7 miles away from the sub- 
station. Inbound signals were measured at an intermediate point and 
at the substation. Computer models were made of these tests for 
comparison. 

BACKGROUND 

Carrier signals were transmitted at a capacitor bank location shown 
on Figure 1. Inbound signals were monitored at Carlton Rd., and the 
substation. These points are shown in Figure 1. Signals were trans- 
mitted by using the power factor correction capacitors already in place. 

A diagram of this transmitting arrangement is shown in Figure 2. At 
Carlton Rd., individual phase voltages were monitored. Figure 3 shows 
a diagram of this voltage monitoring arrangement. 

It was not possible to readily measure individual phase voltages 
at the substation during the field measurements, although the individual 
phase currents were available. Therefore, the phase currents were used 
for model verification purposes at the substation. Figure 4 shows a 
diagram of the measurement apparatus at the substation. 
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RESULTS 


Measurement of inbound signal propagation was made under 
three conditions: 

1) No coupling and measurement apparatus at Carlton Rd. 

2) Coupling and measurement equipment attached at Carlton 
Rd. No tuning inductors attached. 

3) Coupling and measurement equipment attached at Carlton 
Rd. Tuning inductors attached. 

Therefore, under condition (1), no voltage measurements were 
obtained at Carlton Rd. Figures 5, 6, and 7 show comparisons of 
measured and computed values under condition (1). The ratio of the 
individual phase currents at the substation to the transmitted 
voltage at the substation is shown. (Transadmittance) 

Figures 8, 9, and 10 show the ratio of the individual phase 
currents at the substation to the transmitted voltage under condition 
2. Also, under this condition, the individual phase voltages were 
monitored at Carlton Rd. The ratio of the individual phase voltages 
to the transmit voltage at the capacitor bank is shown In Figures 11, 
12, and 13, respectively. 

Under condition 3, the tuning inductors were attached to the 
coupling network. Figures 14, 15, and 16 show the ratio of the 
Individual phase currents. at the substation to the transmitted 
voltage at the capacitor bank. Figures 17, 18, and 19 show the 
measured and computed voltage transfer ratios for the individual 
phase voltages between the capacitor bank and Carlton Rd. 
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FIGURE 19 





CONCLUSIONS 


When making calculations of inbound signal propagation, it was 
necessary to find some representation for the admittance looking back 
Into the substation itself. The substation admittance matrix was 
obtained by field measurements. These measurements were not performed 
at the same time the inbound propagation measurements were made. Since 
the 34556 feeder shares a common substation bus with other distribution 
feeders, the input admittance of the substation could be quite variable 
due to changing load levels on the other feeders sharing that bus. 
Measurement values are directly affected by the substation admittance 
and have an impact on model verification accuracy. 
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I . Overvi ew 


The sensitivity of PLC propagation on feeders to parameter varia- 
tions was investigated using the distribution network computer model . This 
investigation was concerned with two overall goals. The first goal of this 
study was simply to determine the relative significance of various para- 
meters on the characteristics of the network as well as on the propagation 
of a communication signal. The second and more practice-oriented goal of 
the study was to determine if "set*up" time for computer descriptions of 
actual distribution networks could be significantly reduced. The under- 
lying concept of this goal was that if certain parameters were found to 
be insignificant, variations in them could be ignored, thus reducing the 
complexity of the computer description. 

The study was carried out by repetitive analysis of a generic 
distribution network with each execution of the analysis programs having 
one or more parameters varied. Comparison routines were added to the soft- 
ware so certain values could be examined with respect to the corresponding 
values of a reference network. Differences in the compared values were 
recorded as a percent of the reference value. This representation was 
selected so that the degree of sensitivity to a given parameter variation 
could quickly be judged. Other advantages of this representation include 
the ability to directly compare the sensitivity to variations in para- 
meters that are not related, as well as the ability to compare numerous 
variations of the same parameter. 

The following sections deal with the specifics of the study. 

Section II is a discussion of the values compared by the study to deter- 
mine sensitivity - what values were compared, why were these values chosen. 
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and how were the comparisons made. Section III examines the parameters 
that were varied in the course of the study, detailing why each variation 
was made. Section IV contains comments on the more significant results of 
the study and Section V drawn conclusions from these results. Additional 
information pertinent to this study is included in Section VI. An Appendix 
is attached which includes a description of each network that was studied. 

No in-depth discussion of theory is attempted in this report. It 
is assumed that the reader is somewhat familiar with the engineering prin- 
ciples behind the phenomena associated with the operation of a Power Line 
Carrier (PLC) communication network and that further explanation would 
prove redundant. Focus is made, instead, on results of the study and 
conclusions drawn from these results. Naturally, the aspects investigated 
during the study were influenced to a great extent by their relevance to 
the overall goals as well as by a host of other factors including time 
constraints and applicability to the general case. Above all, it must be 
understood that this is solely a generic demonstration intended to point 
out basic trends in sensitivity. In no way are the results to be con- 
sidered typical of all distribution networks. 

II . Comparisons 

Four properties of the distribution network were selected to form 
the basis of comparison for all parameter sensitivity variations. Two of 
these four, the characteristic admittance (Y ) and the eigenvalues, can be 
thought of as "internal" to the feeder (i.e., involving changes in the E 
and H fields in and around the wires) and thus, independent of variations 
in parameters "external" to the feeder such as transformer loading. The 
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other two properties, driving point admittance and voltage transfer, 
incorporate "internal" as well as "external" parameter variations. 

During evaluation of variations in the so-called "external" parameters, 
it was not necessary to make comparisons of the characteristic admittance 
or of the eigenvalues. The other two properties were, however, compared 
for each and every parameter variation made during the study. 

Each of the four network properties is dealt with separately below. 
The discussions present reasoning for selecting that property, as well as 
a brief description of how the actual comparisons took place. Because a 
direct comparison is relatively straightforward, this report omits lengthy 
details of the comparison routines. 

1. Characteristic Admittance - This complex square matrix is 
dependent only on the line configuration and was chosen 
for comparison for several reasons. First, it strongly 
influences the admittance levels associated with an actual 
line and its external loading. Secondly, via its off- 
diagonal terms, it gives an indication of the degree of 
cross coupling between phases. Finally, being dependent 
only on the line configuration, it, along with the complex 
vector of eigenvalues which are dependent only on the line 
configuration and can therefore be computed independently 
of any network, affords some basis of judging the relative 
significance of different line types. 

The element by element comparison took place in the 
following manner. A square complex matrix was initialized 
with the calculated elements of the characteristic admit- 
tance of the reference line configuration. The elements 
along the matrix diagonal are usually dominated by the 
complex characteristic admittances between each phase and 
the neutral, but also contain some cross coupling 
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admittances. The off-diagonal elements indicate only the 
cross coupling effects between the individual phases. This 
matrix was saved for subsequent comparison. Analysis of net- 
works containing parameter variations produced characteris- 
tic admittance matrices which could be directly compared, 
element by element, to the saved reference matrix. Because 
all matrix elements were expressed in rectangular form 
(i.e., a + jb), the difference matrix, 

difference = Y reference - Y varied 
ooo 

was obtained by direct complex matrix subtraction. Each 
complex element of the difference matrix was divided by the 
magnitude of the corresponding element of the reference 
matrix and then multiplied by 100% to express all differences 
as a percent of the reference value. Also presented, for 
quick reference, are the largest diagonal and the largest off 
diagonal difference percentages. 

See Figure (1) for an example of output from the 
characteristic admittance matrix comparison routine. 

2. Eigenvalues - This complex vector is also dependent only on 
the line configuration. It is of interest since each element 
represents the propagation constant associated with the res- 
pective mode. As such, the real part represents the modal 
attenuation (per meter) and the imaginary part determines 
the mode propagation velocity. It is chosen for comparison 
for obvious and similar reasons as cited previously for the 
characteristic admittance. 

The comparison of eigenvalues was carried out in a 
manner analogous to that used for the characteristic admit- 
tance matrix. Differences in the two routines are simply 
logical progressions from the fact that no matrices are in- 
volved, only complex values. Also, since there were no 
diagonal and off-diagonal elements, only the single largest 
difference magnitude was calculated. An example of output 
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LINE TYPE 103 


.* FREQUENCY CODE 10 


INPUT - CHARACTERISTIC ADMITTANCE MATRIX (Y()> FOR LINE TYPE 103 


2.3327974E-03 -6.6306631 E-04 -3. 797 1 1 59E-04 
3.6818 999E-Q5 -3.6365634E-08 7 . 4458229E-06 


-6„ 63C6635E-04 
-3.636.779 4 E-03 


2.4953793E-03 -6..5S343 I 0E-04 
3.6658397E-05 -6.5C5997 2E-03 


3.7371 159E-04 -6.5884303E-04 
7. 4455 229E-06 -6. 5058055E-08 


2.39 28889E- 03 
3.6784427E-05 


REFERENCE - CHARACTERISTIC ADMITTANCE MATRIX (YO) FOR LINE TYPE 101 
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-7.5721 202E-04 

-4.21 43699E-04 

b6 16V43E-06 

6.644 1 I00E-06 

1 . 16753 23E-05 

572 120! E-04 

2.7149607E-03 

-7.5250554E-04 

6441 103E-C6 

2. 3539 32 7E- 05 

6.5767893E-C6 

2I437C4E-C4 

-7.5250551 E-04 

2.586C260E-03 

I675324E-05 

6.5T:H900E-06 

2.5481 140E- 05 

; r£R£NCS MATRIX (ALL FIGUR 

ft 

ES IN PER CENT) 

- -7.45 

1 2.43 

10. 13 

0.44 

-0.88 

-1 .00 

12.43 

-3.09 

12.45 

-0.83 

0.49 

-0.33 

10.13 

12.45 

-7.47 

-l.CO 

-0.38 

0.44 
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MAXIMUM ON DIAGONAL DIFFERENCE ' 8.10 % 


MAXIMUM OFF DIAGONAL DIFFERENCE 12.48 % 

. • 

FIGURE (1) 
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from the eigenvalue comparison routine may be u.»en in 
Figure (2). 

3. Driving Point Admittance - Another complex square matrix, 
this was chosen for comparison because it provides a complete 
electrical representation of a feeder or branch of a feeder 

as seen from the point of input. As was mentioned before, the 
driving point admittance includes all loading effects on. the 
network, those inherent to the feeder as well as those due to 
external devices or conditions. Accordingly, the entire 
feeder or branch of a feeder may be replaced mathematically 
by this single matrix. 

The driving point admittance is similar in most respects 
to the characteristic admittance, the major exception being 
that it encompasses the total range of loading on the network. 
As with the characteristic admittance, the elements on the 
matrix diagonal represent the individual phase to neutral 
admittances, plus phase to phase cross-coupling admittances, 
while the off-diagonal elements consist of only cross-coupling 
admittances. The two matrices are identical in form, so the 
same comparison routine was used for both. Figure (3) is an 
example of output from the driving point admittance matrix 
comparison routine. 

4. Voltage Transfer - This aspect of network performance was 
chosen for comparison for obvious reasons. Because the trans- 
fer of voltage from the input point to the output point on 
the feeder network is the bottom line measure of performance*, 
this property was used as the major factor in determining 

the sensitivity of the computer model to parameter varia- 
tions. 


♦Because the distribution transformer is most easily viewed in terms of 
its own voltage transfer ratio. 
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INPUT - EIGENVALUES FOR LINE TYPE 103 


6.131021 7E-C6 1 .8990793E-06 2. 1 98635 I E-06 

1 .5964033E-04 I . 789740 1 E-04 1 . S026546E-04 


REFERENCE -EIGENVALUES FOR LINE TYPE 101 


5.7937637E-06 7.9675423E-07 9 . 3040934E-07 

1.3665690E-04 I .72397 1 OE-04 1 . 7264795E-04 


DIFFERENCE VECTOR (ALL FIGURES IN PER CENT) 


0.13 0.64 0.73 

r.49 3.81 4.41 


MAXIMUM DIFFERENCE ELEMENT 4.47 % 


FIGURE (2) 
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LINE TYPE 101 


FREQUENCY C0D= 10 


NETN0 n K ID # II. 


INPUT - D°I V I Nj POINT ADMITTANCE «f*T3TX FOP NET VQO< * || 


I-.6335 I 4QP-03 
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5.36953 45—04 

1 .601 6656 --O 3 
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HR I VINO POI IT ADM I 'IT \\ f C- ’* a TR l '< FOP* N-TNOE< * 10 


1 .65.37 155E-03 -9. 7036034^-04 
I .5169 /34E-0 3 5. 401 334 7—04 


-9. ?5 i H730--0 A 
3.45416 133-0 t 


-3. 79 3605 IE-04 
5.4() |*3344E-04 


1 .656 1 334E-0 7 
I .5956I44E-03 


-9.6731 9 75 E— 04 
5.37101 75=-0 4 


-9. 25 10 722E-04 -3 .6791 3372—04 l .60 7739 72—0 3 
3.4541 6692-04 5.3/191 6 3 E-04 - 1 . 500 3 3392-0 3 


DIFFERENCE MATRIX 

(ALL FIGURES IN 

PER CENT) 

-0.23 - „ 

-1.03 

-1.03 

0.53 

-0.09 

-0. 4 7 

• 

-1.03 

-0.0.9 

- 1 .05 

-0.09 

• 

0.61 

- 0 . 0 ; 

-1.08 

'-1 .05 

- 0 . 1 0 

-0,47 

-0.04 

0.66 

MAXIMUM ON DIAGONAL 

njFFE9EN.CE 

0 » 67 


MAXIMUM OFF DIAGONAL DIFFERENCE 1.13 


FIGURE (3) 
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The voltage transfer was established by applying an input 
of 1.0+ jO.O volts on each of the three phase conductors at 
the root section of each network analyzed in the study. 

Assuming there was no other excitation anywhere along the line, 
the resulting voltage level was calculated at the terminal 
section of each network. The output voltages of the three 
phases were presented as the magnitudes of the output voltage 
phasors. The comparison routine received these values for 
the reference network and compared them to corresponding values 
determined for each altered parameter network. Any differences 
between the output phase voltages and their respective reference 
values were computed as a percentage of the reference magnitude. 
Because the magnitude of the input voltage in every case is 
known, a simple hand calculation can provide the change in 
output voltage magnitude as a percent of the input voltage. 

This property can provide a more stable standard of comparison 
because the percent change is referred to the same value no 
matter which reference network was used. It also gives a 
better feel of how much variation can be expected because of 
a certain parameter variation. For an example of output from 
the output voltage magnitude comparison routine, see Figure (4). 

Also generated during the analysis of voltage transfer were 
plots of the phase voltage magnitudes versus distance from the 
input point. These plots provide a graphic representation of 
standing wave and line attenuation effects on the voltage level 
along the line. Production of the plots was more or less a 
sideline activity performed to enhance the understanding of the 
effect of parameter variations. In many cases they illustrate 
damping due to increased attenuation, or other phenomena such 
as altered standing wave patterns due to changes in line 
length or speed of propagation on the phase wires. A sample 
plot can be seen in Figure (5). In this and all other cases, 
the dotted line is the magnitude of the voltage on phase A, 
the solid line is phase B, and the dashed line is phase C. 
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-i . ro 
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FIGURE (4) 
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III. Parameter Variations 

The parameters that were varied for the purpose of this study 
were chosen in accordance with their relevance to the overall goals. 
Although some limited investigation was made into the compounding of 
effects when two or more parameters are varied, most comparisons were made 
with only a single variation. Doing so maintained a "feel" for sensitivity 
to each parameter that could easily be lost by the compounding of effects. 
In some cases the parameter being dealt with was binary in nature (i.e., 
either "yes" or "no"). In these cases, it was not possible to establish a 
normalized sensitivity. Most variations, however, allowed for simple cal- 
culations to determine the relationship between degree of parameter varia- 
tion and degree of change in the model's computations. 

No discussion of parameter variations would be appropriate without 
first describing the reference used in comparison. In the course of the 
study it was necessary to establish four separate references. The 
reasoning here becomes obvious when the scope of the parameter variations 
is taken into account. First, a reference network was needed for varia- 
tions on art overhead feeder. A second reference was needed when investi- 
gation was made into variations on an overhead feeder approximately one 
quarter wavelength in total length. The great disparity in length between 
the original and the one quarter wavelength networks made comparisons to 
the initial reference somewhat uninformative. Following the same logic, 
two more references were needed when underground cable networks were 
investigated. Although confusing, four reference networks were needed 
to maintain a sense of scale in the sensitivity study. 
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All four reference networks were essentially the same feeder but 
each case was altered to fit the specified application. The basic con- 
figuration common to all four was a synthetic, cascaded line having no 
branches, no three-phase to two- or one-phase conversions, and no varia- 
tions in line configuration anywhere along the feeder. It was broken up 
into 80 equal length sections, each section Iqaded by a 25 KVA distribution 
transformer cyclically connected from phase to neutral. The individual" 
transformer admittances were set at -78.8 db at an angle of -27.3°. The 

line was terminated by a load having equal matrix-diagonal elements of 
-4 

4.77 x 10 mhos at -68.22° and zero off-diagnoal elements. 

The nominal overhead reference network had sections 200 meters long 
for an overall length of 16,000 meters. Each of the 80 sections was loaded 
by one of the distribution transformers. When an approximately one- 
quarter wavelength overhead network was to be investigated, the length 
of each section was reduced to 110 meters for an overall length of 8,800 
meters. This value was established by calculating the wavelength of the 
8.13 kHz PLC signal on the nominal reference network using the computed 
speeds of propagation for each of the three modes. The 8,800 meter length 
was chosen as a round approximation of the three values calculated for one 
quarter wavelength. Both the 16,000 meter overhead and the 8,800 meter 
version were, by unfortunate choice, referred to in the data as Network #10. 

The nominal underground cable reference network was identical to the 
nominal overhead case in most respects. The only differences were the 
obvious change in line configuration and an increase in the temperature 
of the conductors. Whereas the overhead case assumed all conductors to be 
at 25°C, the phase conductors in the underground reference were set at 
75°C with the neutral strands at 45°C. All other network parameters were 
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identical to the overhead version. The length of the one-quarter wave- 
length underground reference network was computed by the same process 
used for the overhead line. It resulted in the network being reduced from 
16,000 to 3,400 meters. The cable reference networks are referred to in 
the data as Networks #9 and #8, for the nominal and quarter wavelength 
versions, respectively. For more information concerning the reference 
networks, see Figure (6) or refer to Appendix 8. 

Now that the reference networks havs been established, the parameter 
variations may be collected into four groups. The groups are divided, in 
most cases, according to the intent of the parameter variation. The first 
group consists of variations in overhead configuration intended to 
establish the relative sensitivity to such changes. The second group 
included all loading variations intended to determine loading effects. 
Investigations into "set up" time reduction constituted the third group 
while cable network variations were separated into the fourth. Obviously 
the groups are far from mutually exclusive since the results of variations 
in; one group may be of significance to the intentions of variations in 
another. Divisions are made along logical lines for the clear presenta- 
tion of results . 

A. Overhead Geometry Variations - The following variations were made in 
the geometry of the overhead feeder: 

1. Phase Conductor Diameter - Three variations were made, both 
increasing and decreasing the diameter from its reference 
value. The investigation was carried out because changes in 
wire diameter are a common occurrence on real-life feeders. 

If such changes proved to be insignificant to a communication 
signal, the set up procedure could be simplified without 
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sacrificing significant model accuracy. If, however, effects 
were notable, the degree of sensitivity would be of interest. 

Analysis was performed on both the nominal system and the 
quarter wavelength version. Perhaps a comment on the quarter 
wavelength investigation is appropriate. The investigation 
was performed to ascertain the sensitivity of quarter wavelength 
networks with relation to longer networks. It was suspected that 
the unique standing wave patterns existing on feeders approxi- 
mately one quarter wavelength long could possibly magnify sensi- 
tivity to a given variation. To test this, the same variation 
was analyzed on the two networks and the resulting sensitivities 
were compared. In short, the quarter wavelength analysis was 
performed to see if sensitivity was magnified by standing waves, 

2. Neutral Conductor Diameter - One variation of the neutral conduc- 
tor diameter was made on both the nominal and the one quarter 
wavelength networks. Sensitivity to the neutral was suspected 

to be somewhat less than sensitivity to the phase conductors 
so less attention was afforded this parameter. Reasoning 
behind this variation was much the same as that behind the 
phase conductor diameter variations. 

3. Phase Conductor to Phase Conductor Spacing - One variation was 
made in which the spacing was increased. Both the nominal and 
the quarter wavelength networks were analyzed. 

4. Phase Conductor to Earth Spacing - One variation was made on 
both network lengths. The neutral conductor to earth spacing 
was held constant so the phase conductor to neutral spacing 
was increased along with the phase conductor to earth spacing. 

5. Neutral Conductor to Earth Spacing - One variation was made in 
which the phase conductor to earth spacing was held constant as 
the neutral to earth spacing was increased. As with all over- 
head geometry variations, analysis was performed on both the 
nominal and the one quarter wavelength networks. 
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6. Miscellaneous Conductors - Three variations were made involving 
miscellaneous conductors being added to the pole below the 
phase conductors. First one and then two additional neutral 
conductors were added, all assumed to be at ground potential. 
Finally, two open circuited secondary conductors were simulated. 
These were also assumed to be essentially at ground potential. 

7. Triplex Neutral Conductor - Analysis was performed to determine 
the effect of changing the neutral conductor from simple ACSR 
to Triplex cable. 

8. Spacer Cable Configuration - One investigation was made into the 
sensitivity of the model to a radical change in overhead geometry. 
The standard cross arm configuration common to all other analysis 
was altered to a so-called spacer cable or Hendrix cable arrange- 
ment. See Figure (8-1) in Appendix 8 for an illustration of this 
geometry. 

B. Loading Variations - Investigation was made into the sensitivity of 
the model to variations in loading on the network. These variations 
involved distribution transformer loading, the diagonal load at the 
end of the network, and a series of combinations to provide informa- 
tion on the compounding effects of multiple variations. It was known 
from previous work that loading would have a pronounced effect on the 
output of the model so this group of variations was performed to 
determine the degree of sensitivity. Very little "set up" time 
reduction was foreseen. 

1. Terminal Load - Analysis was carried out on both the nominal 
and one quarter wavelength networks with the diagonal load at 
the end of the line removed. The line was left as an open 
circuit. 
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2. Distribution Transformer Loading - Many variations were per- 
formed involving the loading due to distribution transformers. 
These loads were doubled, cut in half, or completely ignored. 

It should be noted that when these loads were doubled or cut 
in half, it was the load as seen by the primary that was 
being altered, not the physical loads on the secondary side 

of the transformers. Because of the large series leakage re- 
actances that were present, doubling the secondary load would 
not necessarily double the load seen by the primary. All 
references to altering loads in this report should be considered 
in the context of the load on the primary side. Most, but not 
all of these variations were performed on the one quarter wave- 
length network as well as the nominal network. 

3. Terminal Load and Transformer Loading Combination - One investi- 
gation was made into the combined effects of removing the terminal 
load and doubling the distribution transformer loading. The 
degree to which one variation would counteract the other was of 
interest. The variation was made on both network lengths. 

4. Transformer Loading and Overhead Geometry Combination - In order 
to investigate the effects of changes both "internal" and 
"external" to the feeder, the entire series of overhead geometry 
variations was repeated with the distribution transformer loading 
cut in half. It was considered possible that geometry variations 
might be more significant if the damping effect of the transformer 
loading was reduced. It was also of interest to note if the 
variations compounded one another. This series of variations 

was conducted on the nominal length network only. 

C. "Set-up" Time Reduction Variations - The largest number of variations 
were those concerned mostly with the streamlining of the network des- 
cription procedure. This time-consuming and exhaustive procedure 
requires the detailing of every physical aspect suspected to effect 
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the propagation of a PLC signal. Individual investigation was made 
into a fair number of these parameters in an attempt to determine 
if excessive attention was being paid to insignificant parameters. 

It was hoped that in many instances, prudent engineering assumptions 
could replace costly precision determinations without a significant 
loss of model accuracy. 

In addition to variations in specific parameters, changes in 
"set up" procedures were also looked into. Specifically, the proce- 
dure of creating a new section at each distribution transformer, even 
though there is no change in the feeder. Also investigated was the 
compounding of effects when "set up" variations interacted with 
loading variations. Many of the "set up" variations were analyzed on 
the one quarter wavelength network in addition to the nominal length 
version. 

1. Aggregation - Two variations were made involving the reduction of 
the total number of sections in a network description. Both were 
demonstrations of aggregation of several shorter sections into a 
single longer one. Because all sections in the synthetic test 
network were of the same configuration, it was theoretically 
possible to create one long section with the corresponding 
included distribution transformers located at the end of this 
long section. Variations made in the study, however, were 
restricted to the more reasonable possibilities of reducing the 
number of sections to one half and one third of the original 80. 

The aggregation process was accomplished conveniently in 
the model to calculate the effect by doubling (or tripling) of 
the length of every other (or third) section and reducing the 
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in-between sections to zero length.* Because 80 is not evenly 
divisible by three, it was necessary to merely double the length 
of the next to last section in the second case. In both cases, 
the same overall length was maintained, as was the number and 
loading of the distribution transformers. Although in actuality, 
the number of sections in the network description was not reduced, 
the effective number of sections was indeed reduced and the 
technique of aggregation was demonstrated. The same technique 
was tested on a real life feeder description. The results of 
that test are discussed later in this report. Aggregation 
analysis was performed on both nominal and quarter wavelength 
networks . 

2. Temperature - A series of five temperature variations was made on 
both length networks. The temperatures investigated covered a 
wide but possibly low range of conceivable conditions. This study 
was made to determine if any great concern need be placed on the 
value of the conductor temperature. 

3. Earth Resistivity - Four variations were made covering a range of 
values likely to be encountered in most areas of the country. 

Both quarter wavelength and nominal cases were investigated. 

4. Ungrounded Neutral - All investigations performed other than this 
one assumed the neutral wire to be essentially at ground poten- 
tial. This analysis removed that assumption and was done on both 
nominal and one quarter wavelength networks. 

5. Aggregated Networks with Transformer Admittances Doubled - The 
aggregation study described earlier was repeated with the loading 
of the transformers on the feeder doubled. The intent was to see 
if variations were greatly compounded. 


*In actual practice, a simple modification to the model and input file 
procedure would be made to make this aggregation more convenient. 
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6. Total Length - The length of each section of the nominal net- 
work was reduced to 190 meters for an overall reduction of 5 % 
in the length of the network. To obtain a 5 % increase in total 
length, the length of each section was upped to 210 meters. 

It was hoped to establish the sensitivity of the model to the 
accuracy of line length measurements. 

In a sense, the same variation can be thought of as a varia- 
tion in the signal frequency rather than in total line length. 

If the line loading is held constant, the output voltage for a 
given input will depend on how many wavelengths are contained 
on the line at steady state. A given line can be considered to 
be N wavelengths long, and as long as N remains constant, it 
doesn't matter whether the line length varies slightly or the 
wavelength varies slightly. In this case, the speed of propa- 
gation was constant so a change in frequency would result in a 
change in wavelength. Mathematically, 

£ = NX 

X = S P/f 

£ = NS P oh NS = £f 

~T P 

where £ is the length of the line, Sp is the speed of propa- 
gation, X is the wavelength, and f is the signal frequency. 

Since Sp is a constant and the object is to hold N constant, 
an increase in £ must be accompanied by a decrease in f and 
vice versa. From this it can be seen that increasing the 
line length while holding the signal frequency constant is 
essentially the same as holding the length constant and 
decreasing the signal frequency, if the same number of wave- 
lengths is to be contained. 

D. Cable Network Variations - Many of the variations made on overhead 
feeders were repeated for underground cables. Although the total 
number of variations did not approach that of the overhead networks. 
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those performed were comprehensive enough to indicate cable sensi- 
tivity to most parameters of interest, It should be noted that 
every variation made on the nominal underground network was also 
analysed on a one quarter wavelength underground network. As with 
the overhead case* this was done to determine any possible sensiti- 
vity magnification on networks of that length. 

The parameter variations performed on the underground feeders 
were analysed for much the same reasons as the corresponding overhead 
variations. The cable variations are considered separately from the 
overhead feeders because of the large differences in the characteris- 
tics of the two configurations. The very low characteristic Impedance 
of the underground cable as compared to overhead lines dramatically 
alters the sensitivity to variations, Also, the lowered speed of 
propagation radically changes the standing wave patterns on feeders 
of the same length, Separate consideration simplifies the presenta- 
tion of results and conclusions made from those results, 

1, Earth Resistivity - Similar to the overhead case, the earth 
resistivity was varied over a reasonable range of values, 

Temperature - Cable feeder descriptions allowed for separate 
specification of phase conductor and neutral strand tempera- 
ture, In all three variations, the neutral strands were set 
30° lower than the phase conductors. The temperatures 
evaluated were quite high when compared to those analysed for 
the overhead case. This reflects the thermal insulation provi- 
ded by the cable l s electrical Insulation and by the ground. 
Analysis at overhead temperatures would not have reflected 
conditions on existing cable feeders. 


3. Aggregation - The analysis of simplified network descriptions 
was handled as it was for overhead networks * 

4. Transformer Loading - Again* this was handled as it was for 
overhead lines. Admittances were doubled, cut In half, and 
neglected altogether, 

IV, Results 

The results of the parameter sensitivity study are most clearly 
presented when the four major groups of variations are further subdivided 
according to the parameter varied, Within each of these subdivisions, 
discussion Is made concerning the effects on the various network pro- 
perties compared for the study. Two of the four properties* the charac- 
teristic admittance and the driving point admittance, are matrix quanti- 
ties and thus present values both on and off the matrix diagonal . When 
the method of signalling employed involves placing identical signals on 
all three phases, it Is the elements on the matrix diagonal that are of 
greatest concern. These elements include phase to neutral admittances 
as well as cross-coupling admittances whereas off-diagonal elements 
contain only cross-coupling terms. In general, the phase to neutral 
admittances are an order of magnitude or so greater than the cross- 
coupling terms and tend to dominate, 

This section attempts to make a capsulated presentation of signi- 
ficant results. In order to prevent the inundation of the reader by 
detailing thousands of percentages and related minutla, every effort 
has been made to reduce measures of sensitivity to as few numerics as 
possible. In each case, the percent differences in the magnitudes of 
the three phase output voltages were averaged to provide a single value 
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to be used for comparison, For other properties, the largest difference 
magnitude (on-diagonal for matrices) was used because it represents 
the limiting case. Whenever possible, a normalized sensitivity is 
presented. This value is defined as the percent difference in the 
output voltage magnitude divided by the percent difference in the 
parameter. Both the numerator and the denominator are expressed as a 
percent of their reference values. Mathematically, the normalized 
sensitivity equals 


AV 

“V 


out 

out 


x 100% 


y- X 100 % 


where 


V ou t = reference value of output voltage magnitude 

AV out = change in value of the output voltage magnitude 

P = reference value of the parameter 

AP = change in value of the parameter 

Because changes in the output voltage magnitudes are by far 
the most significant indicators of parameter sensitivity, extra 
attention was given to this property. In addition to figuring the 
change in output voltage magnitude as a percent of the reference 
output voltage magnitude, a calculation was performed to determine 
the change in output voltage magnitude as a percent of the input 
voltage magnitude. In each case analyzed, the input voltage was 
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1 .0 /0° volts, so. 


- /V out' - /V out!/ 


in 


'V 


x 100% 


where V out 
v 

v out 


/V / - /V 7 
- xl00i 

* -'''out/ - /V out’' x 10W 
= reference output voltage 
= output voltage with a parameter varied. 


There was also a normalized sensitivity calculated from this 
property. In mathematical terms, 

/V .ut' - /v .ut‘' x 10OT 
££ x 100* 

where all values are defined as in earlier equations. 

Again, let it be noted that these results were obtained on 
highly idealized networks. Since the intent is merely to indicate 
general trends, the results will be presented in as general terms as 
possible. More specific relationships may be determined via direct 
analysis of the data. 

A. Overhead Geometry Variations 

1. Phase Conductor Diameter - The three variations in phase 
conductor diameter consisted of 21 .9% and 38.0% decreases 
and a 22.5% increase. Results were tabulated for all four 
of the properties compared by the study - characteristic 
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admittance, eigenvalues, driving point admittance, and 
voltage transfer. 

Results indicated that sensitivity to phase conductor 
diameter was not consistent across the range evaluated. The 
two decreases in diameter produced normalized sensitivities 
referred to the input voltage of 0.119 and 0.108, values 
similar enough to be considered constant when computational 
errors are included. On the other hand, the lone increase in 
phase conductor diameter indicated a sensitivity of only 0.047, 
far less than half the values determined in the other variations. 
The sensitivities of the other three properties pretty much 
followed the pattern established by the output voltage magni- 
tude. In every case the sensitivity of the property was 
significantly less for the phase conductor diameter increase 
than for the decreases. See Figure (7) for the specific 
numerical values. 

One quarter wavelength investigation revealed no in- 
creased sensitivity to phase conductor diameter variations. 
Normalized sensitivities referred to the input voltage assumed 
nearly the same three values determined in the nominal evalua- 
tion 0.117, 0.108, and 0.50, respectively. The characteristic 
admittance and the eigenvalues were not recomputed for the 
quarter wavelength study because neither was length dependent. 
Exact results of the one quarter wavelength study may be seen 
in Figure (8). 

2. Neutral Conductor Diameter - The only variation made was a 37.2% 
reduction in diameter which resulted in a normalized sensitivity 
with respect to the input voltage of just 0.018. Changes in 
every category were slight even though the parameter variation 
was quite large. Just as with the phase conductor diameter, 
no change in sensitivity was noted on the quarter wavelength 
system. Here the final sensitivity measured was just 0.017, 
virtually the same as on the nominal network. See Figures (7), 
and (8) for specific results. 
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3* Phase Conductor to Phase Conductor Spacing - A 10.01 increase 
in this parameter produced a normalized sensitivity referred 
to the input voltage of only 0.067, One quarter wavelength 
evaluation resulted in the similar value of 0,071. As with 
all overhead geometry valuations, the characteristic admit- 
tance and the eigenvalues were computed only once even though 
both nominal and one quarter wavelength networks were analyzed. 
Numerical results for the two cases may be seen in Figures (7) 
and (S). 

4. Phase Conductor to Earth Spacing - This parameter was Increased 
10,0$ while the neutral conductor to ground spacing was main- 
tained. (In effect, the phase conductor to neutral conductor 
spacing was increased approximately 52.51). This resulted in a 
final sensitivity measure of -0,355, The negative indicates that 
an increase in spacing results in a decrease in output voltage 
magnitude. (This value is reduced to approximately -0.067 if 
the phase to neutral spacing is considered to be the parameter 
varied). Again, there was no particular increase in sensitivity 
when the one quarter wavelength network was analyzed. Normalized 
sensitivity referred to the' input voltage for that case was -0.360, 
Refer to Figures (7) and (8) for other specific ratios and per- 
cent changes. 

5. Neutral Conductor to Earth Spacing - Effectively the opposite of 
Increasing phase conductor to earth spacing, a 10.01 increase was 
made while phase to earth spacing was held constant. (In effect, 
the phase conductor to neutral conductor spacing was decreased 
approximately 42.51). The result was a normalized sensitivity 
referred to the input voltage of 0,427 (or alternately 0.10). 

As was standard for overhead geometry variations, no Increase 
in sensitivity was noted when the one-quarter wavelength net- 
work was analyzed. Results on that network were 0.443 (or 
(0.104). See Figures (7) and (8). 
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Miscellaneous Conductors - Three variations were made on both 
nominal and quarter wavelength networks that involved the 
addition of miscellaneous conductors to the overhead configu- 
ration. Because these variations are binary in nature, no 
normalized sensitivities could be computed. Instead, compari- 
sons will be based on the change in output voltage magnitude, 
referred to either the reference output value or the input 
value. 

The introduction of a single miscellaneous conductor at 
ground potential resulted in a 12.7% increase in the output 
voltage level, as compared to the reference level. Because 
the reference output voltage magnitude was approximately half 
of the input magnitude, the increase, as compared to the input 
voltage, was approximately half of that compared to the output 
voltage, specifically 5.75%. This condition was a general trend 
with all overhead geometry variations. The other miscellaneous 
conductor variations involved two conductors being added to the 
pole. These resulted in voltage increases of 15.2% and 14.5%, 
referred to the reference output voltage magnitude. The same 
increases were 6.87% and 6.56% when referred to the input 
voltage magnitude. 

One quarter wavelength analysis produced much the same 
results as the longer network. Comparisons here should be made 
with the change in output voltage magnitude as a percent of the 
input voltage. This is due to the fact that the reference out- 
put voltages are different on the nominal and the one quarter 
wavelength networks, thus, changes as a percent of these values 
should not be directly compared across networks. Comparisons 
within each length network are legitimate. More complete pre- 
sentation of results may be seen in Figures (7) and (8). 

Triplex Neutral Conductor - Changing the neutral conductor to a 
triplex cable was another binary type variation and again, no 
normalized sensitivity could be obtained. Comparing output 
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FIGURE (7) 



voltage magnitudes showed a 5.31% increase over the reference 
value, or a 2.40% increase as compared to the input voltage. 

The same values determined on the one-quarter wavelength network 
were 4.50% and 2.29%, the latter of which should be used for 
comparison and is within 5% of the value for the nominal network. 
As a general note, all neutral variations produced less of a 
voltage drop as a percent of input on the quarter wavelength 
network than on the nominal length network. All of the quarter 
wavelength values were, however, within 5% of the nominal 
values, so although consistantly less, they were not siqnificantly 
less. See Figures (7) and (8) for details. 

8. Spacer Cable - This change in configuration was so radical it 
was impossible to establish a degree of change. Accordingly, 
no normalized sensitivity was determined. As was expected, 
there was a significant change in the output voltage, up 21.4% 
on the nominal and 22.7% on the quarter wave network. Both 
figures are increases as a percent of the input voltage magni- 
tude. Changes in the other properties were comparably great. 

See Figures (7) and (8) for the specific values. 

Loading Variations 

1. Terminal Load - The equal element diagonal load at the end 
section of the feeder was removed on both the nomimal length 
and one quarter wavelength networks. Because this was another 
binary type of variation, no normalized sensitivity could be 
obtained. Comparing the change in output voltage as a percent 
of input voltage, shows a significant difference between the 
nominal and quarter wave cases. The nominal network showed 
a 9.98% increase while the quarter wave network jumped 17,1%. 

No measurements were made of characteristic admittance or 
eigenvalues because this variation would not effect those pro- 
perties. See Figure (9) for more detailed results. 
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Distribution Transformer Loading - The results of variations in 
this parameter established a rather significant link between 
transformer loading and voltage level on the feeder. The doubling 
of the transformer loading as seen by the primary side can be 
thought of as a 100% increase in loading. Similarly, halving 
the loading can be considered a 50.0% decrease. These varia- 
tions can thus result in normalized sensitivities. Ignoring 
the loading altogether is a binary variation and cannot be 
expressed in normalized terms. Transformer loading was doubled 
and neglected on both nominal and one quarter wavelength net- 
works. Halving was performed on only the nominal length network. 

Doubling of the transformer loading resulted in a normalized 
sensitivity with respect to the input voltage of -0.297 for 
the nominal length and -0.282 for the quarter wavelength case. 

As before, the negative value indicates that loading increases 
result in voltage decreases. On the other hand, decreasing the 
loading by 50% resulted in a normalized sensitivity of -0.716, 
indicating a sensitivity twice as great as the earlier results. 
This carries through to the other properties to the same general 
degree. As with analysis of the terminal load variation, 
characteristic admittance and eigenvalues were not compared 
because they would not be effected by this variation. 

The analysis performed with all transformer loading ignored 
produced very large differences in the output voltage. More 
importantly, the voltage increase on the one quarter wavelength 
network was many times larger than that on the nominal network. 

The change in output voltage magnitude as a percent of the input 
voltage was 268% for the quarter wave network and 59.3% for the 
16,000 meter version. Because removal of distribution trans- 
former loading eliminates the vast majority of damping effects, 
the measured output voltage is almost purely dependent on the 
input voltage and the particular standing wave pattern on the 
line. As expected, an undamped one quarter wavelength pattern 
greatly magnified the output voltage level. More complete 


3-32 


results for all distribution transformer loading variations 
may be found in Figure (9). 

3. Terminal Load and Distribution Transformer Loading Combination - 
Both the quarter wavelength and the nominal length network 
results indicated that combining these two variations allowed 
them to cancel each other to some degree. Using AV out /V.j n x 100% 
as the yardstick for comparison, nominal length analysis resulted 
in a 25.7% decrease in output voltage. This compares to a 29.7% 
decrease with the transformer loading doubled and a 9.98% increase 
with the terminal load removed. Although it is not a linear 
combination of the components, there is a definite interference 
between the two. The same pattern was repeated for the quarter 
wavelength results but the interference had a slighter effect. 
Results were a 23.3% drop in output voltage as compared to a 17.1% 
increase and a 28.2% decrease when the terminal load and trans- 
former loading are dealt with separately. A smaller drop would 
have been expected due to the larger increase from removing the 
terminal load, assuming the two components counteracted linearly. 
Obviously, they do not. Figure (9) supplies additional results 
that may be of interest. 

4. Transformer Loading and Overhead Geometry Combination - The 
results of this series of variations continued along the lines 
of the previous combination. Because there was more than one 
variation being made, it was not possible to establish a normali- 
zed sensitivity that would have had any worthwhile meaning. 
Instead, the change in output voltage magnitude as a percent of 
input voltage magnitude was the standard of measure. In every 
case, this ratio indicated a significant degree of interference 
between the effects of the two variations. As an example, 
consider the case in which the distribution transformer loading 
was halved and the phase conductor diameter was decreased 21 .9%. 
These variations considered separately produced output voltage 
changes of +35.8% and -2.60%, respectively. Taken in tandem. 
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FIGURE (9) 









result was a 32.7% increase in voltage, clearly indicating some 
sort of additive effect. Details of the results for this series 
may be found in Figures (9) and (10). Values for the characteris- 
tic admittance and eigenvalues reflect the fact that overhead 
geometry variations altered these properties. The change in 
driving point admittance reflects both the altered characteris- 
-tic admittance and the reduced admittances of the transformer 
loading. 

C. "Set-up” Time Reduction Variations 

1. Aggregation - If the number of sections in the network is con- 
sidered to be the parameter varied in these cases, a normalized 
sensitivity can be obtained for the aggregation analysis. The 
two variations were thus a 50.0% reduction and a 66.3% reduction 
in the effective number of sections. Results indicated that 
normalized sensitivities with respect to the input voltage magni- 
tude were 0.009 and 0.014 for the two respective cases. Quarter 
wavelength analysis produced the similarly insensitive values 

of 0.011 and 0.015. Characteristic admittance and eigenvalues 
were not compared for these variations because they would not 
be effected by the number of sections in the network. See 
Figures (11) and (12) for more complete results. 

2. Temperature - The sensitivity of the model to variations in the 
temperature of the overhead phase conductors was rather firmly 
established. The results of all five temperature variations 
indicated a normalized sensitivity with respect to the input 
voltage of only -.001. Repetition of these same variations on 
a one quarter wavelength system produced values too small to be 
meaningfully recorded. Similarly, the driving point admittance 
was largely unaffected by temperature. The characteristic 
admittance and the eigenvalues were not recorded even though 
they would have been altered by temperature variations. Aside 
from the fact that obtaining comparisons would have been quite 
difficult for various reasons, the sensitivities computed indi- 
cate that differences would have been far too small to have been 
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FIGURE (10) 








of great interest anyway. Actual results tv iu ^ comparisons 
may be found in Figures (11) and (12). 

3. Earth Resistivity - The four earth resistivity variations per- 
formed on the quarter wavelength and the nominal length networks 
resulted in a nearly constant value for the normalized sensiti- 
vity referred to the input voltage magnitude. All eight values 
were within the range of -0.01 to -0.02. All other measures 
showed similarly low sensitivity. As with the temperature varia- 
tions, comparisons of characteristic admittance and eigenvalues 
were not carried out even though they would have been effected 

by changes in earth resistivity. The same reasons can be cited. 
See Figures (11), (12), and (13) for the results of comparisons 
made . 

4. Ungrounded Neutral - Analysis performed with it no longer assumed 
that the neutral was everywhere at ground potential resulted in 
very slight changes in the output voltage magnitudes. Because 
this was a binary variation, normalized sensitivities could not 
be computed. In the absence of this, the change in output 
voltage magnitude as a percent of input voltage magnitude was 
used as a yardstick. Results showed a .297% increase on the 
nominal length network and a .423$ increase on the quarter wave- 
length network. In either case, the increase was small. See 
Figures (12) and (13). 

5 . Agg r egati on Networks With Distribution Transformer Admittances 
Doubled - Because two variations were made at once, no norma- 
lized sensitivity was calculated. Using the change in output 
voltage as a percent of input voltage as a comparison, the 
results of these variations could pretty much have been pre- 
dicted by looking at the results of the component variations when 
considered separately. Doubling transformer admittances de- 
creased voltage by approximately 29$. Aggregation reduced output 
voltage approximately 1$. When combined, voltage reductions were 
30.0% and 30.2% for networks with total sections reduced 50% and 
66.3%, respectively. The figures were 28.5% and 28.8% reductions 
for one quarter wavelength analysis. These are also in line with 
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the values found from quarter wavelength analysis of the single 
parameter variations. Again, a slight interference effect can 
be noted as the two reductions combine. The doubled loading, 
of course, is the much more significant of the two variations. 
Complete results may be found in Figures (11) and (12). 

6. Total Length - Increasing and decreasing the overall length of 
the nominal network resulted, surprisingly enough, in voltage 
drops for both cases. This is indicative of an alteration in 
the standing wave pattern rather than a mere change in the point 
of measurement along that pattern. Study of the voltage vs. 
distance plots supports this. At any rate, the normalized sensi 
tivities referred to the input voltage were .718 for the 5% 
decrease and -.988 for the 5 % increase. These values show an 
almost one-to-one ratio of A V to AP and point out that total 
length is a parameter of significance. See Figure (13) for 
additional results. 

Cable Network Variations 

1. Earth Resistivity - The sensitivity of underground cable net- 
works to variations in earth resistivity was even less than 
that of overhead feeders. Normalized sensitivity with respect 
to the input voltage was, at most, a mere -0.0015. This figure 
includes both nominal length and one quarter wavelength results. 
Typically the results were values too small to meaningfully 
record the actual digits. The results, as presented in Figures 
(14) and (15), indicate that for all intents and purposes, under 
ground cables are impervious to variations in earth resistivity. 

2. Temperature - The results presented here were calculated consi- 
dering only the change in temperature of the phase conductors. 
Because the feeder descriptions of underground cables allow for 
separate specification of phase conductor and neutral strand 
temperatures, ths percentage differences of the two were not 
equal even though the number of degrees difference from the 
reference values were equal. When the percentage difference 
was needed in a calculation, the value for the phase conductors 

was used. It should also be noted that percentage differences 
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FIGURE (11) 
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were determined in degrees Centigrade. The alternate procedure 
of determining differences in degrees Kelvin would have resulted 
in much smaller percentage differences for the same variations 
and numerical accuracy would have been lost. 

Results showed a marked increase in sensitivity as compared 
to overhead lines. The sensitivity measure was only -0.001 for 
the overhead case but jumped to an average value of -0.094 for 
the nominal length cable and -0.036 for the quarter wavelength 
case. Also notable is that the one quarter wavelength sensiti- 
vity is over three times the nominal length sensitivity when 
both are computed referred to the input voltage. When sensiti- 
vity is normalized considering the change in voltage as a per- 
cent of the reference output voltage, the values indicate that 
the nominal network is the more sensitive. This phenomenon 
is due to the voltage magnification caused by the standing wave 
pattern on the quarter wavelength network. The output voltage 
on that network was in the vicinity of three times the value of 
the input. This greater than unity voltage transfer and the 
associated magnification of sensitivity are examples of the 
behavior that was the object of the quarter wave investigation. 
Complete results of the temperature variations on both length 
networks may be obtained from Figures (14) and (15). 

3. Aggregation - Results of these variations again pointed to 
quarter wavelength magnification. Normalized sensitivity as 
referred to the input voltage was infinitesimal on the nominal 
length network. Exact values were 0.0003 for a 50% reduction 
in the number of sections and 0.0004 for a 66.3% reduction. 

The same variations on a one quarter wavelength cable resulted 
in values of 0.0067 and 0.01, respectively. Although still 
insignificant to the operation of a PLC network, they none- 
theless indicate the magnifying effect of the standing wave 
pattern. See Figures (14) and (15) for additional details. 
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4. Distribution Transformer Loading - Unlike the overhead lines, 
variations in this category had only a small effect on the 
voltage transfer of the modeled networks. As with the other 
cable variations, there is evidence of tremendous sensitivity 
magnification on the one quarter wavelength network. The sensi 
tivities with respect to the input voltage magnitude calculated 
on the nominal length network were slight but constant. A 100% 
increase in loading resulted in value of -0.0134 while a 50% 
decrease in loading gave a sensitivity of -0.0137. Of course, 
these values reflect a lot of averaging and round-off and are 
not really accurate to that many decimal places. Even though, 

a high degree of consistency is indicated. The corresponding 
values from the quarter wavelength analysis were -0.254 and 
-0.288. Most notable is the way the quarter wave magnification 
has made the loading a parameter of some significance to commu- 
nications on this network. 

Completely ignoring the transformer loading could not be 
expressed in terms of normalized sensitivities. Looking at 
AV out^in x again points up the magnification of sensiti- 
vity due to a total length of essentially one quarter wave- 
length. A relatively small increase of only 1.37% was exper- 
ienced on the nominal length network when all transformer 
loading was removed. This was multiplied to 29.3% on the 
quarter wavelength feeder. Again, the parameter variation 
has become of some significance to communication. More 
complete details may be obtained from Figures (14) and (15). 

5. Total Length - The sensitivity determined by these variations 
indicated that length was relatively important to the model 
predictions. The sensitivity calculations produced values of 
-0.228 for a 5% decrease in line length and 0.358 for a 5% 
increase. These values establish total length as the most 
significant parameter on the nominal length underground cable. 
See Figure (16) for complete results. 
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FIGURE (15) 



FIGURE (16) 








V. Conclusions 


Conclusions drawn from the results of the parameter sensitivity 
study are presented in two groups. These groups correspond to the two 
overall goals of the study - to determine the sensitivity of the model 
to various parameters and to evaluate possibilities of reducing "set up" 
time. Accordingly, the conclusions will be based on two slightly 
different criteria. The first goal was more or less a pure knowledge 
endeavor while the second was specifically practice-oriented for a PLC 
communications network. It would, therefore, be possible for the model 
to be highly sensitive to a parameter, making that parameter of great 
interest for the first goal, yet because the parameter varies only over 
a narrow range, it is insignificant to a PLC communications network. 

Actual instances such as described here occurred during the course of 
the study. 

Due to the fact that performance of a PLC communications network 
would not suffer serious adverse effects unless the signal level dropped 
by several dB, the voltage transfer ratio of the network would necessarily 
have to be altered by 40 or 50 percent before any great significance could 
be assigned to a parameter variation. On the other hand, the sensitivity 
of the model was better judged by the ratio of change in output to change 
in parameter. By this determination, the voltage transfer did not neces- 
sarily have to change by a great amount to establish a high degree of 
sensitivity. These two criteria formed the basis for the conclusions of 
this report. 


3-47 


A. Sensitivity of the Model to Various Parameters 


1. Overhead Geometry - The model showed greatest sensitivity to 
radical changes in the overhead configuration. The change 
to spacer cable from the reference cross-arm caused the 
voltage transfer ratio to jump to nearly 150% of its reference 
value. Approximately the same change was recorded on the 
quarter wavelength system so it can be concluded that there is 
no increase in sensitivity to a change of this nature. From 
the results of this variation, it should be determined that any 
change in the basic overhead format (i.e., cross-arm to armless, 
armless to spacer cable, etc.) will result in a significant 
change in the model predictions if the lengths of such configu- 
rations are significant. 

Also causing a notable increase in the voltage transfer 
ratio is the addition of miscell aneous conductors. Simply 
adding a wire or two bumped the voltage to 110% or 112% of its 
initial value, slightly more for quarter wavelength analysis. 

Even changing to a triplex neutral wire improved the ratio, 
but only by 5% or so. 

Of the variations for which normalized sensitivities could be 
computed, it turned out that changes in the phase to earth and 
neutral to earth spacing resulted in the greatest changes in the 
voltage transfer. If, he /er, these variations are considered to 
be phase conductor to neutral conductor spacing variations, the 
normalized sensitivity value is greatly reduced. The modelled 
network proved to be relatively insensitive to phase to phase 
spacing and neutral conductor diameter. Only phase conductor 
diameter showed signsof slight significance to the model. 
Sensitivity to this parameter was such that small variations 
(less than 30% or so) could easily be ignored but larger ones 
begin to have noticeable effects on the model computations. 

In virtually every case, the model showed little or no increase 
in sensitivity when the network length was reduced to 
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approximately one quarter wavelength of the communication 
signal. This is largely the result of the transformer 
loading, which remained constant through all of these 
variations and was responsible for most of the damping of 
the standing waves via signal attenuation. 

Generally speaking, unless there is a change in the basic 
geometric pattern of the overhead configuration, the model will 
be only slightly sensitive to variations. Changes in the 
neutral wire are mostly insignificant. Miscellaneous conductors 
will improve the voltage transfer ratio by a few percent. Phase 
conductor to neutral conductor spacing and phase conductor dia- 
meter should be watched because sizable variations will have 
effects on the model calculations. 

2. Loading - In contrast to overhead geometry, the model proved to 
be quite sensitive to changes in the loading along the line. In 
every single variation made involving loading, the voltage trans- 
fer ratio computed by the model was significantly altered. In 
the many combined variations involving loading, the effects of 
the loading variation tended to swamp out those of the other 
variation made in the combination. Because the leakage reac- 
tances of the distribution transformers constitute a significant 
portion of their total loading effect, loading variations on the 
secondary side will not necessarily produce a variation of the 
same magnitude when seen from the primary side. This tends to 
dampen variations but nonetheless, loading can vary widely. The 
variations analyzed attempted to cover a broad range by taking 
a reference load and ranging it from half that value to twice 
that value. It would be fair to say that at the limits of this 
range, one could expect the voltage transfer ratio to have 
changed by approximately 30% (at least on the synthetic test 
network). Logically, as the loading increases the voltage trans- 
fer ratio decreases and vice versa. 
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The quarter wavelength analysis performed for loading varia- 
tions provided a wealth of information. Most importantly, it 
pointed out the tremendous damping of the standing wave pattern 
caused by the loading. When the loading, as seen on the primary 
side of each distribution transformer was doubled, virtually the 
same change in the voltage transfer ratio was noted on both net- 
work lengths, When the loading was reduced or eliminated, the 
quarter wavelength network showed signs of greatly Increased 
sensitivity. As the damping of the standing wave pattern 
lessened, the voltage magnification became more and more pro- 
nounced. The extreme case occurred when the distribution trans- 
former loading was ignored altogether. In this case, the voltage 
transfer ratio increased on the quarter wavelength network by 
over four times the amount that it increased on the nominal length 
system. The quarter wavelength analysis demonstrated how much of 
the standing wave suppression is due to transformer loading. To 
a lesser degree, the terminal load on the line also damped the 
standi nq waves. This was also demonstrated by this study. 

The conclusions drawn from study of loading variations ar' 
that loading is the most significant parameter to overhead lines. 
This sensitivity is caused by the damping effect the loading has 
on standing waves. Also, there appears to be a saturation point, 
beyond which model sensitivity to loading increases drops off 
dramatically. This was indicated by the fact that there was no 
significant magnification in sensitivity on the quarter wavelength 
network when the reference loading was doubled. The standing 
wave pattern had reached a point where it was so damped that it 
was no longer responsive to quarter wavelength magnification. 

3. "Set-up 1 * Time Reduction - This group of variations can be broken 
up so the parameters varied can be dealt with separately. These 
parameters were the number of sections in the network, the 
temperature of the conductors, the earth resistivity, and the total 
length of the line. It can safely be said that the model is 
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completely insensitive to all of these parameters except total 
length. Twenty-six variations were conducted with the first 
three parameters and not one indicated any signficant sensiti- 
vity in the model. Also, analysis with the neutral no longer 
assumed to be at ground potential pointed up no sensitivity to 
that parameter. The one quarter wavelength analysis simply 
pointed out that any effects due to these parameter variations 
were insufficient to overcome the suppression of standing waves 
caused by transformer loading. 

Total length variations were of some consequence, however. 
The sensitivity of the model to length variations was due to 
the alteration of the standing wave pattern caused by changing 
the length. The variations were not great enough for a change 
in the sum total of the characteristic impedance to be of any 
Importance. The only logica 1 conclusion is that model sensiti- 
vity was a direct consequence of standing wave pattern sensiti- 
vity. As was discussed earlier, this variation can also be seen 
as a variation in signal frequency. Thus, the propagation will 
be sensitive to some degree to variations in either parameter. 
Since the standing wave pattern will be different for each 
specific length of line, signal frequency, and loading level, 
it would be impossible to state beforehand how much of an effect 
a certain error in line length estimation or signal frequency 
will have. The sensitivity of the propagation will vary with 
each Individual case. 

4. Cable Networks - Because of the very low characteristic impedance 
of underground cables, in most cases they proved to be far less 
sensitive to parameter variations, particularly in the case of 
the nominal length network. On that network, even the neglecting 
of all transformer loading increased the voltage transfer ratio 
by little more than 1%. This near total insensitivity reflects 
the fact that the impedances of the distribution transformers 
and other loads are so high when compared to that of the cable, 
that they are seen almost as open circuits. Variations in these 
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external loads are of little consequence when the model is 
dealing with cables. The greatest sensitivity noted was to 
variations in the temperature of the phase conductors. Even 
this sensitivity was slight. Sensitivity to earth resistivity 
and aggregation was too small to accurately determine. 

One quarter wavelength analysis pointed out that even 
though sensitivity was reduced, it still could become significant 
if standing waves are such that voltages are magnified. Under 
these conditions, the loading on the line produced sizable damping 
effects, as did increases in the temperature. It would be safe 
to conclude that under quarter wavelength magnification effects 
the propagation becomes sensitive to loading and temperature 
variations . 

Variations in the total length of the nominal underground 
cable showed a small but possibly significant degree of sensiti- 
vity. Results showed this sensitivity to be about one-third of 
that on the overhead line version of the simple test network. 

It can be argued along the same line used in the overhead case 
that this sensitivity can vary greatly depending on each specific 
case. Therefore, it is difficult to make conclusions on its 
importance based on the data available. 

B . Possibilities of Reducing "Set-up" Time 

Almost every variation made had some implication with regard to 
"set-up" time reduction. The most direct of these were, of course, the 
aggregation demonstrations. Many other variations indicated that the 
model did not necessarily require a highly accurate value for certain 
parameters - a single, rough estimate for the entire network would 
provide accurate enough results. Obviously, reducing the amount of 
accurate information needed will reduce the complexity of the "set up" 
procedure and thus the time and cost. For this reason, conclusions will 
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be drawn from results of many variations not directly intended to 
reduce "set-up" time. 

Here we encounter the question - How accurate is accurate 
enough? As was discussed earlier, a change in voltage transfer of 
perhaps 40% is the minimum variation significant to a PLC communications 
network. Only a handful of variations made during the study had such 
an impact. The combination variations indicated that there was an 
additive effect when two variations occured together. It is this effect 
that must be accounted for when the standards of accuracy are established. 

1. Overhead Geometry - These variations indicate that changes in 
the neutral conductor can largely be ignored. Sensitivity of 
the model to these variations was so small that a single, 
standard neutral may be assumed with little loss of accuracy. 

Only two aspects of the neutral need be considered. If the phase 
to neutral spacing varies by a large amount, or if additional 
neutral wires are present, then these variations should be 
accounted for. This is true only if the variation is present 
on a significant portion of the network. If it exists on only 
100 meters of a 10,000 meter line, it could easily be ignored. 
However, 3,000 meters out of 10,000 might be significant enough 
for consideration. 

Along similar lines, only significant lengths of Tine con- 
taining miscellaneous conductors or phase conductor diameter 
changes need be bothered with. Actually, small changes in con- 
ductor diameter can be ignored completely no matter how much of 
the line they effect. Even very short lengths of radical geometry 
changes can safely be passed over. Any change that exists for 
more than 5 % of the total length should be considered significant, 
however. Variations in phase to phase spacing can safely be 
ignored. 
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2. Loading - A sizable variation in loading on the line may be 
enough to significantly effect the voltage transfer all by itself. 
Because of the sensitivity here, it may be well to conduct a 
systematic loading variation study over a plausible range of 
relative loading to determine a worst case for each modeled 
feeder. 

This is particularly true on quarter wavelength networks 
because of the greatly enhanced sensitivity to loading, 
especially light loading. It- would be safe to say that no 
variation in loading greater than 10% or 15% should be ignored. 

3. "Set-Up 11 Time Reduction - The most promising results of these 
variations indicated that aggregation of several shorter sections 
into one long one was an effective way of simplifying network 
descriptions without introducing any significant loss of accuracy. 
As was alluded to earlier, the demonstration of this technique 

was taken one step further. Analysis was performed on an actual 
computer description of a real life network. Every possible 
aggregation was made and the new description was reanalized with 
results compared to the original. The original network may be 
found in the data as network #72. The aggregated version is 
network #73. Results of this real life test case showed that the 
number of sections in the network was reduced over 36.8%, yet 
the voltage magnitudes at the terminal section were changed by an 
average of only 0.187% of the input value. Clearly a significant 
simplification can be achieved with no great less of model 
accuracy. It is important to note that total length and trans- 
former loading were held exactly constant during these demon- 
strations. 

The other results of this group of variations also indicate 
that descriptions can be further simplified. Overhead networks 
are for all intents and purposes completely insensitive to 
variations in conductor temperature and earth resistivity. This 
allows basically one default value to be used for almost every 
network. Results will be sufficiently accurate and there need 
be no undue concern over the exact values for temperature or 
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earth resistivity. Also, results point out that model calcula- 
tions are not very dependent on the assumption that the neutral 
is at ground potential everywhere. It can be inferred that the 
model description will be accurate even if this assumption is 
totally untrue. 

The only variation that showed a necessity for accuracy 
in description was the total length of the network. Because the 
standing wave pattern is so easily upset, it is somewhat impor- 
tant that overall length be judged within\+- 5% of its actual 
val ue. 

4. Cable Networks - Conclusions regarding cable networks are not so 
concrete. As was pointed out before, variations that indicated no 
sensitivity at all on the nominal length network resulted in 
sizable sensitivity values when analyzed on a quarter wavelength 
line. Consequently, a knowledge of how long the network is in 
terms of signal wavelengths might be helpful in the process of 
deciding which variations can be safely ignored. No matter how 
long the feeder is, the total length is a parameter that should 
be documented within _+_5% or thereabouts. As with overhead 
lines, the length of the entire network is somewhat important and 
must be accurately determined if the model is to produce valid 
results. 

The loading on the line and the temperature of the phase 
conductors are parameters that varied in significance depending 
on whether or not analysis was performed on a quarter wavelength 
line. At that length, they were significant enough to warrant 
attention. Any variation of 20% or more that could be determined 
should be considered significant. At longer lengths, these para- 
meters fade in significance and no tremendous attention need be 
paid to them. Specification of a default temperature value would 
perhaps be accurate enough, and any realistic loading value would 
suffi ce . 
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Finally, aggregation and earth resistivity had almost no 
effect on the model. This leads to the conclusions that aggre- 
gation can be carried out successfully on underground cables, 
and that a default value for earth resistivity is accurate 
enough for any computations. In the general case, "set up" 
procedures for underground cable networks can be streamlined 
even more than for overhead lines because modal sensitivity is 
generally lower. 

By far the most important conclusion reached by this study is that 
aggregation is a valid method for reducing network description complexity 
without adversely effecting results. In conjunction with other simpli- 
fications, this will significantly reduce the time and effort now being 
expended for accurately detailed descriptions. By determining variations 
in certain parameters to be of little significance, the number of different 
line types and configurations can now be greatly reduced. Earth resis- 
tivity and conductor temperatures need not be accurately measured. These 
and other simplifications will go a long way toward making model descrip- 
tions easier to deal with and yet just as accurate. 

Establishment of sensitivities for the model also aids in actual 
operation of the programs by enhancing the user's feel for the tool he 
is utilizing. Although these results were obtained on a highly idealized 
network, they shed light on other networks. Basic trends can be identi- 
fied and applied elsewhere. Unusual phenomena may be understood when 
knowledge of the model's sensitivities is available. AH in all, appli- 
cation of the model may provide results that can now be interpreted in 
a more meaningful way. 


VI. Comments on Numerical Round-off Error 


As with any process involving a great deal of computation, some 
error is introduced by rounding-off or truncating numbers. This was 
initially thought to be significant in some instances, so an evaluation 
was made. This evaluation involved the analysis of a 10,000 meter, 3- 
phase line, having a single 25 kva transformer load. Analysis was per- 
formed both as a line with 100 sections and as a line aggregated into 
one section. The following results indicate the overall voltage transfer 
matrix computation was subject to, in most cases, <_ 1* error due to 
round-off. This indicates that only a very small part of all error is 
due to numerical round-off error. For all practical purposes, the error 
is of no great significance. See Figures (17) and (18) for the results 
of this experiment. 
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Attachment to DITNAP Memo #26 (Cont'd) 


Voltage Transfer Matrix Obtained by Matrix-Matrix Multiplication 
100 Section Line 


2.2205505E-00 

-1.0050040E-01 

T.5582932E-01 

-7.4422400E-02 

1 .6295392E-01 
-8.661 561 4E-02 


1 .2949408E-01 
-6.1201580E-02 

2.1916808E-01 

-8.8912366E-02 

1 .3002984E-01 
-6.8101263E-02 


1.9955678E-01 
-1 .5897823E-01 

2.0303388E-01 
-1 .6863907E-01 

2.3463748E-00 
-3.0640301 E-01 


Voltage Transfer Matrix Single Section Line 


2.2207349E-0Q 
-1 .01 66037E-01 

1.5608705E-01 

-7.5870808E-02 

1 .6379893E-01 
-8.9224265E-02 


1 .2963878E-01 
-6.2126620E-02 

2.1918909E-00 
-9.0067551 E-02 

1.3070214E-01 
-7.1059561 E-02 


1 .9850288E-01 
-1 .5408064E-01 

2.0209737E-01 
-1 .641 9739E-01 

2.3349881 E-00 
-3.0308890E-01 


Output Voltage for a Dirac-Delta Voltage Drive on Phase 1 
100 Section Using Matrix-Vector Multiplication 


2.22073E-00 1 .56087E-01 1 .63794E-01 

-1.01660E-01 -7.58708E-02 -8.92234E-02 


100 Section Using Matrix-Matrix Multiplication 

2.2205505E-00 1 .5582932E-01 1 .6295392E-01 

-1.005004E-01 -7.4422400E-02 -8.5515614E-02 

Single Section Line 

2.220735E-00 1.560870E-01 1.637989E-01 

*1 .016604E-01 -7.587081 E-02 -8.922426E-02 


FIGURE (17) 
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Attachment to DIFNAP Memo #26 
Input Admittance 100 Section Line 


Input Admittance 100 Section Line 
DB-Angl e 


-46.23 

-58.62 

-65.86 

88.58 

-87.66 

-77.21 

-58.62 

-45.90 

-58.72 

-87.66 

88.83 

-86.04 

-65.86 

-58.72 

-45.59 

-77.21 

-86.04 

83.52 


1 .21 28379E-04 4 .7941084E-05 

4.8796294E-03 -1 .1714650E-03 


Rectangular 
1 .127641 7E-04 
-4.96801 88E-04 


4.79411 34E-05 
-1 .1 714569E-03 

1 . 1 2761 63E-04 
-4 . 9681 725E-04 


1 .0354070E-04 
5.0669665E-03 

8.0014310E-05 

-1.1565352E-03 


8.0014984E-05 
-1. 1565391 E-03 

5.9263794E-04 
5.2191951 E-03 


Input Admittance 1 Section Line 


Input Admittance Single Section Line 
OB-angle 


-46.23 

-58.62 

-65.86 

88.58 

87.66 

77.21 

-58.62 

-45.90 

-58.72 

-87.66 

88.83 

-86.04 

-66.86 

-58.72 

-45.59 

-77.21 

-86.04 

83.52 


Rectangular 

1.21 2841 8E-04 4.79411 97E-05 1 .1276341 E-04 

4.8796358E-03 -1 .1714630E-03 -4.9680814E-04 


4.79411 95E-05 
-1 . 1 71 4629E-03 

1.1276339E-04 

-4.9680826E-04 


1.0354079E-04 

5.0669705E-03 

8.001 4676E-05 
-1 .1565380E-03 


8-001 4670E-05 
-1 .1565380E-03 

5.9263896E-04 
5.2191 989E-03 


FIGURE (18) 
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TASK 4 

TRANSFORMER MODELING 


R. C. Rustay 
J. T. Gajjar 
R. W. Rankin 
A. M. Dunham 


FOREWORD 

This chapter reports on the activities associated with modeling 
of two winding transformers in the context of use at PLC frequencies. 
Distribution transformers are included in this category by assuming 
them electrically symmetric and having balanced secondary loads with 
the net result they become effectively two winding transformers. 

PURPOSE 

The purpose behind the effort to obtain a suitable computational 
model was to obtain an automated software procedure (not requiring 
manual intervention) for predicting two-way propagation through the 
transformers, i.e., from primary to secondary and secondary to primary. 
Also, with such a model, the effect of different levels of secondary 
loading (deterministic and random, i.e., MONTE CARLO) could be easily 
accomplished. From a system engineering point of view, PLC propagation 
through transformers represents a critical item in system performance 
and having a physical modal is very useful for understanding these limi- 
tations and the effects of possible fixes such as tuning and bypass 
procedures . 

Transformers of interest are: 

Distribution Transformers (assuming electrical symmetry 
and balanced secondary loading) 
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Single Phase Ratio Bank Transformers 

Three Phase Transformer Banks (Three Single Phase 
Transformers) including open delta banks. 

Single Core Three Phase Transformers were not considered. 

Synopsis of Phase I Transformer Modeling Activities 

Transformer modeling was initiated in Phase I of this project. 

At that time, little previous apriori knowledge existed regarding the 
behavior of distribution transformers at PLC frequencies and, in 
particular, when the various higher order levels of dlatributed/stray 
capacity became important. Also, it was originally intended that 
models be developed for frequencies to 500 kHz. Therefore, it did not 
appear that an analog type lumped parameter model could be developed 
and an associated parameter data base established within the time, 
available manpower and funding constraints associated with the Phase I 
effort. Therefore, it was decided to base modeling on measured "Y" 
parameters which would be frequency dependent (and of course KVA and 
other constructional features) but otherwise completely general. This 
approach was encouraged by several references on the subject*. 

Therefore, concurrent activities of measurement, and model 
development and software implementation were initated. An extensive 
set of "Y" parameter measurements** were made on 5, 10, 15, 25, 50 KVA 
distribution transformers for frequencies to 500 kHz. Subsequently, 
when attempting to use these measured "Y" parameters, significant 
inconsistencies were discovered. After much deliberation, it was 


*$ee pgs. 635-641 "Magnetic Circuits and Transformers", MIT EE Staff 
John Wiley & Sons. 

**Six "Y" parameters plus an open circuit impedance. 


hypothesized that the "¥*' parameters must be measured/determined with 
an accuracy and precision that was not realized and would probably not 
be practical even with more sophisticated procedures. 

By the time this conclusion was reached, only limited resources 
and time remained. In order to temporarily overcome this problem, a 
procedure was implemented using a file containing measured primary 
driving point admittances (for specified level of secondary loading). 
This allowed predictions to be made on the feeder for the specified 
loading implied in the transformer file being used. Prediction "down 
to" the secondary would be a second manual step using again measured 
voltage transfer characteristics. This "two step" procedure was not 
in any way deficient as an approximation but did require manual inter- 
action, and did not admit random "MONTE CARLO" procedures. 

Synopsis of "Inter Phase" Transformer Modeling Activities 

During the approximately one year interval between the conclu- 
sion of PHASE I contract effort and the commencement of this PHASE II 
contract, the General Electric Co. sponsored extensive response type 
measurements on distribution transformers, single phase power trans- 
formers, and various single core three phase transformers. These 
measurements were made at Corporate Research and Development with the 
cooperation of the Niagara Mohawk Power Company [NMPC] who provided 
access to the various transformers, and Mobile Radio Department, both 
General Electric operations. The measurements included bi-directional 
driving point admittance and voltage transfer ratios as a function of 
frequency (limited range) and various loadings (primary and secondary). 
It was from these measurements that adequacy of a single RLC lumped 


parameter model for frequencies up to, say, 50-100 kHz, was suggested. 

Synopsis of PHASE II Transformer Modeling Activities 

Still faced with need and solution for an adequate transformer 
model, this activity was continued as Task 4 in the PHASE II contract. 
Four modeling approaches were considered: 

1) All "Y" parameter model (as in PHASE I) with greater 
measurement precision. 

2) "Y" parameter derived algebraic model 

3) Table look-up of measured responses 

4) Simplified lumped parameter model 

Approach 1) was rejected because of the reasons cited previously and 
the uncertainty whether sufficient measurement precision could be 
obtained. Approach 3) was rejected (after 2 and 4 appeared feasible) 
because of the large number ,)f interpolating dimensions (frequency, 
magnitude of loading, primary and secondary, phase of loading) and the 
prohibitively large amount of measurements required. 

During the above deliberations, an algebraic basis was developed, 
corresponding preliminary measurements were made, and a computer program 
written to establish the feasibility of approach 2). Appendix 3 contains 
analyses relevant to this approach. The key difference between this 
approach and the "Y" parameter approach 1) was to assume that: 

a) at frequencies of interest, distribution transformers are 
electrically symmetric (two secondary windings), 

b) a balanced secondary load (no basis exists for assuming 
otherwise for purposes of PLC propagation modeling), 

and to replace some of the "Y" parameters with other more readily and 

algebraically less sensitive measurements such as open circuit voltage 
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transfer ratios. As a result of laboratory measurements on several 
transformers, a set of derived parameters was selected which were 
relatively insensitive to transformer construction. The resulting 
algebraic model was implemented in a temporary test (computer) program 
which very successfully predicted responses corresponding to measure- 
ments. The feasibility of this approach appears to be proven. Like 
"Y" parameters, these derived parameters would be frequency and KVA 
dependent, but independent of loading, and would predict bi-directional 
driving point admittances and voltage transfer ratios. It does not 
appear that this model has any frequency limitations (other than 
measuring them at such frequencies). It does not, of course, represent 
a physical model . 

Measurements (see last page of Appendix 3) made during the inter 
phase interval strongly suggested that a simple RLC lumped parameter 
model might be adequate for frequencies limited to, say, less than 50- 
100 kHz. A test computer program was written to evaluate the model 
and quite successfully duplicated the measurements cited above. On 
this basis, subject to the limitations on upper frequency, the simpli- 
fied lumped parameter model also became a viable contender. It should 
be noted that it is realized a more complicated lumped parameter model 
could be derived to extend this frequency limitation, but doing so at 
this time, is not considered necessary. 

Measurement Technique for Phase II 

In determining the performance of a transformer at PLC frequencies, 
four characteristics are significant for the purpose of this program and 
these include; voltage transfer ratio from primary to secondary [^ S /Vp], 


and from secondary to primary [ v P/Vs]» along with the admittance 
looking into the primary [Ypri] and the admittance looking into the 
secondary [Y ] from line to line. All of these quantities must be 
known as a function of primary and secondary loading and the frequencies 
of interest. 

Several distribution transformers were brought into the laboratory 
from the NMPC facility at Seneca St. in Schenectady, NY. A test set-up 
was developed that provided the desired data with sufficient accuracy 
and in an expeditious manner. Figure 1 is a block diagram of the test 
set-up showing the equipment configuration used in all of the measure- 
ments. Terminals 1 and 2 always provide the transformer excitation 
and terminals 3 and 4 are used to monitor either voltage or current. 

The Quantech Wave Analyzer is used solely to provide a ramp 
function, and this ramp function linearly varies the frequency of the 
Lin/Log Sweep Generator from 3 to 100 kHz and, in addition, is used 
to drive the x-input of the x-y plotter. The period of the ramp is 
a front panel adjustment on the Quantech and is set for 500 sec. so 
as to assure adequate response time for the measurement equipment. 

The attenuator shown on the VCG input of the Lin/Log Sweep Generator 
adjusts the frequency range of the sweep generator. 

The output of the sweep generator drives the Bogen power ampli- 
fier which has an impedance matching transformer on its output. Thus, 
the output can be either floating as when driving the secondary line- 
to-line, or have one side tied to ground as when driving the primary. 

The Hewlett-Packard Gain Phase Meter is used to measure the 
magnitude and phase of a voltage-to-voltage or a current- to-voltage 
ratio where the A input is always the reference and is always a 
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Test Set-Up for Transformer Measurements 
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voltage. The response time of the analog outputs of the gain phase 
meter is set to 0.02 sec. by putting the frequency range switch in 
the 100 Hz-lmHz position. There is no digitizing involved in genera- 
ting the analog outputs. 

The 2 analog outputs of the gain phase meter, magnitude and phase, 
are used as inputs to the two y-drives of the x-y plotter. Because of 
the skew on the two pens of the plotter, there is a frequency error 
between the magnitude and phase of about 500 Hz [1 small division on 
the graph paper with the 100 kHz range]. That is, when a magnitude 
value is read at a certain frequency, the corresponding phase value 
should be read at a frequency which is 500 Hz less. 

Figure 2 shows the connections to terminals 1-4 of Fig. 1 which 
are used to measure the magnitude and phase of the voltage transfer 

Vr 

ratio from primary to secondary, S /Vp, of the subject transformer. 

Figure 3 shows the connections used in obtaining Vp /V s . 

Figure 4 shows the connections used in determining the admittance 
looking into the primary, Y . , as a function of secondary load, R g . 

Figure 5 shows the connections for measuring the admittance 
looking into the secondary, Y sgc , as a function of the primary load, Rp. 

The magnitude tolerance of the ratio (B/A) configuration is 
given as +1.5 db for the gain phase meter with input amplitudes of 
greater than -64 dbv [Rj 0.6 mv] and + 2.5 dbv with input amplitudes of 
less than -64 dbv. Similarly, the phase angle error is not more than 
10° in our frequency range with input amplitudes of greater than -44 dbv 
[% 6 mv}. It was judged that these tolerances were in keeping with the 
scope of this program. 
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♦Tektronix Probe ** Tektronix 7603 Oscilloscope with 

Model 6054A 10X 7A12 Dual Trace Amplifier 

[Used as a Differential Amplifier] 
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Transformer Construction 

1 ■■ - ■ - J - 1 ■ 

Two types of transformer construction are in common use, and 
these are classified as core type and shell type. When the magnetic 
circuit takes the form of a single ring encircled by two or more groups 
of primary and secondary windings distributed around the periphery of 
the ring, we have a core type transformer. A simple core-type transformer 
is shown below. 



The term "core type" implies that the coils are cylindrical and 
concentric with an outer winding over an inner winding. 

When the primary and secondary windings take the form of a common 
ring which is encircled by two or more rings of magnetic material distri- 
buted around its periphery, we have a shell type transformer. A simple 
shell -type transformer is shown below. 



The term "shell -type" denotes large pancake coils which are 
stacked or inter-leaved to make primary-secondary groups. 

The characteristic features of the core- type transformer are a 
long mean length of magnetic circuit and a short mean length of windings 
and vice versa for the shell -type. Apparently the choice of which 
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construction to choose is largely a matter of manufacturing facilities 
and individual preference. 

According to the distribution transformer people at General 
Electric, in the 10-100 kva distribution transformer range. General 
Electric Company and RTE make core types and all other domestic manu- 
facturers make shell types. 

Both types were measured during the course of this program because 
it was initially felt that there might be some significant difference in 
their performance at PLC frequencies . No such difference was observed. 

Transformer Symmetry 

The assumption has been made in the transformer modeling program 
that the two halves of the secondary are symmetrical [and have balanced 
loads] so that we can, in effect, deal with a two winding transformer. 

Measurements were made to justify that assumption. The results 
of these measurements are shown in Figs. 6-9. Figure 6 shows the voltage 
transfer ratio, Vp /V s , of a 25 kva General Electric transformer (core- 
type], when the secondary is driven LI -N in one case and L2-N in the 
other case. This measurement is made essentially as shown in Fig. 3, 
with'Rp = 11a, except that the differential amplifier of Fig. 3 is not 
needed here. Figure 6 shows that the magnitude of ^P/V s at a frequency 
of 5 kHz differs by about 4.5 db. Beyond 50 kHz, the two magnitudes 
tend to converge. Because of the way the measurement was made, the 
phase angle associated with ^P/V s is 180° different in each case. The 
phase angle difference at 5 kHz (neglecting the 180°) is about 3% and 
at 50 kHz is about 6°. The maximum angular difference occurs at about 
65 kHz and is on the order of 24°. 
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Figure 7 shows ^P/V s for the 25 kva General Electric transformer 
again, but this time comparing the case of driving 1/2 of the secondary 
as opposed to driving the full secondary. It can be seen that, up to 
about 60 kHz, the difference both in magnitude and phase is trivial. 

Figure 8 shows the y 12 parameter of a 25 kva Colt transformer 
(shell type) with the secondary current measured in the LI -N leg and 
the L2-N leg. The y 12 parameter is defined as the transfer admittance 
looking into the primary with the secondary shorted and given by the 
formula: 

y 12 * 7p— [sec shorted] 
pri 

In Figure 8, I was measured in the Ll-N leg, in the L2-N 
leg, and in the L1-L2 leg. Again, in the frequency range up to 50 kHz, 
the maximum difference in magnitude of y^ 2 did not exceed about 3 db 
and the angular difference did not exceed about 3°. No attempt will be 
made to explain the performance of y^ at frequencies above 60 kHz. 

Figure 9 shows the behavior of y^ of a 25 kva Colt transformer 
when each of the primary bushings [2 bushing transformer] are connected 
to the secondary center tap in turn. The term y^ is defined as the 
input admittance looking into the primary with the secondary short 
circuited and is given by the equation 

I_ ri 

y.. * [sec shorted] 

v pri 

It can be seen that the point of minimum magnitude of y^ has 
changed frequency by about 2500 Hz with a corresponding change in the 

* 

point where the phase angle goes through zero degrees. 
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Finally, the errors incurred as a result of assuming symmetry 
in the two halves of the secondary of a distribution transformer are 
of the same order as those of the instrumentation previously discussed, 
and again, were judged to be in keeping with tha scope of this program. 

Response Measurements 

Response measurements were made on the 25 kva Colt transformer 
and the 25 kva GE transformer, referred to as the laboratory transformer; 
and the results of these are shown in Figs. 10-17. These curves repre- 
sent the kind of information the computer model must be capable of 
generating for each of the transformers associated with the feeder 
section to be modeled. These measurements were made of unenergized, 
as regards 60 Hz power, transformers and with simulated resistive 
primary and secondary loading. 

In regard to secondary loading, most of the response measurements 
were made with the distribution transformers at about 1/2 load. The 
value of load resistor (R^) for full load is determined by the formula 

V 2 

R. _ v sec(L-L) 

L ' kva 

2 

which, for a 25 kva transformer, becomes \ ~ 25 ,000va = 2,2 ohms * 

Primary loading of a distribution transformer in the PLC fre- 
quency range is a function of many things, including type of construc- 
tion, i.e., overhead or underground, number of other transformers on 
the primary, the secondary Toads of other transformers on the primary, 
power factor convection capacitor banks, etc. Our experience has shown 
that the magnitude of the primary impedance varies from about 10 ohms 
[for some underground distribution] to several hundred ohms 
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[for some overhead distribution]. 

Referring to Fig. 10, this shows the magnitude and phase of the 

V s 

voltage transfer ratio from primary to secondary, 77=- . It can be seen 

V v p 

that at about 5 kHz, |7p| is approximately -33 db which is an attenua- 

P 

tion nearly equal to the transformer turns ratio of 31.75 

[20 log 3I775 = -33.04 db]. With increasing frequency |p-| decreases 

P 

and in the frequency range from about 20 kHz to 60 kHz, this decrease 

v s 

is around 6 db/octave. The phase relation for Tp is nearly 0° [or 180° 

, P 

depending on tl>e measurement probe placement] at the lowest frequency 
and then becomes about -66° at 50 kHz, indicating significant reactive 
components . 

v s 

The rp curves for the 25 kva GE transformer, Fig. 14, are almost 
P 

identical to those for the Colt just described. 

The Y . curves of Fig. 11 for the 25 kva Colt transformer exhibit 

pri 

a parallel resonance phenomenon at about 13.5 kHz. The minimum admit- 
tance [maximum impedance] that this 1/2 loaded transformer presents to 
the line is -82 db or 79ymhos. This corresponds to 12.6 kohms of 
impedance. The minimum impedance of this transformer is 1.26 kohms 
and occurs at the maximum frequency (100 kHz). At 100 kHz, Y . is 
almost completely capacitive and has been since about 50 kHz, as 
shown by the phase angle curve. A curve of capacitive susceptance has 
been superimposed on Fig. 11 to show that Y ^ becomes increasingly 
capacitive above 50 kHz and to show the general shape of capacitive 
susceptance curves when plotted on a linear frequency scale. The 
value of capacitance used to plot this curve was C = 1.28 nf. This 
capacitance susceptance curve ca;i be used as an overlay template as 
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an alternative procedure to estimate the equivalent lumped capacitor. 

Aga-in, the Y pri curves of Fig. 15 for the 25 kva 6E transformer 

are very similar to those for the Colt* 

V 

Figure 12 shows the ^ characteristic of the Colt transformer 

s 

for a primary load resistance of 10 and 100 ohms. It can be seen that 

V 

the magnitude curves of the ratio decrease at 6 db/octave up to 

s 

about 50 kHz and further, that the differences between the magnitude 

curves for the 10 ohm and 100 ohm values of primary resistance is 

exactly 20 db over the entire frequency range. It will be shown later 

that the leakage reactance of the transformer essentially forms a 

constant current source when compared to the value of R p transformed 

to the secondary *nd that therefore is proportional to R over a 

w s p 

wide range of R„. 

p V 

The phase curves of indicate that V p and V $ are about 90° 
out of phase up to 50 kHz and in phase [or 180° out of phase] at 
about 100 kHz. 

V 

Again, the curves of the 25 kva 6E transformer. Fig. 16, are 
v s 

very similar to those for the Colt. 

Figure 13 shows the Y curves for the 25 kva Colt transformer 

w Cw 

for a variety of primary resistance [R p ] values ranging from 0 to 

infinity. There is minimal difference in the magnitude and phase 

curves when R„ varies from 0 ohms to 1000 ohms. When R = 0, there 
P P 

is a series resonance at 15 kHz. For values of R of 10,000 and 
1000 ohms, the magnitude of Y decreases at a 6 db/octave rate over 

JW V 

the entire frequency range and the phase angle curve shows Y $ec to be 
mainly inductive over the frequency range. Inductive susceptance 
plotted on a linear frequency scale has the same form as an inverted 
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capacitive susceptance. 

Figure 17 shows t 
GE transformer and the Colt transformer are again similar. 


Figure 17 shows that, for R p * 10 ohms, the Y sec curves of the 


Algebraic Model 

From Appendix 3, the four equations of interest, assuming symme- 
trical secondary windings, line to line secondary measurements and 
drive, and balanced secondary loads are: 


s _ 


= 2 


2 

912 


_2 

y 12 


-1 


where Y 2 = admittance of load on each half of secondary 

v s 

g ^2 = Tpwith secondary open-circuited 
P 

I, 

y 12 * w=-with secondary short circuited 
P 

I n i 1 Y ? 

Y ori = f = t~ + C y n - r~ ] [ — h— ] 

pn V M »/ - y 12 


9 12 


where -J— = with secondary open-circuited 


Z 11 V p 


y ll = v W1th second ary short-circuited 
P 


!e. 

v s 


^21 


1 + Y 


1 

y n 


where Y-j = admittance of primary load 
V 

9 21 = v with P rimary open-circuited 
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The quantities y-] 2 » J — » yy| and 9 2 i were measured, as a 
function of frequency, for both the 25 kva Colt and GE transformers. 
Figures 18-27 show these characteristics. 

V V 

For the 25 kva GE transformer, the value of y-, Y . , and 
- P s 

Y sec were computed at several frequencies using the above four equations 

with data from Figs. 23-27. The results of these computations have 

been plotted on the measured characteristics of Figs. 14-17. 

It can be seen that using the algebraic model with measured 
values of open and short circuited parameters provides reasonable 
correlation with actual measured performance. To use the algebraic 
model in feeder modeling work would require the measurements and 
storage of g-^* y-| 2 > y^ anc * ^21 ^ or transformers of interest 
(at least some representative types) and at all frequencies of interest. 
This is a large amount of data to obtain and to store. 

The advantage of this type of model is that it does not depend 
on a valid lumped parameter RLC model and therefore could be used at 
any frequency providing the necessary measurements were available. 


Lumped Parameter Model 

As mentioned in Appendix 3, earlier work by several parties had 
indicated that a simple "RLC plus ideal transformer" model could 
adequately represent distribution transformers in the frequency range 
of 5-50 kHz, again assuming electrically symmetrical secondaries and 
balanced secondary loads. 





























Figure 21 









Figure 22 
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From this earlier work, the analysis presented in Appendix 3 
and a study of the response measurements of Figs. 10-17, the model 
shown below was developed. 



where N * turns ratio [31 .75] 

V = L-L secondary voltage 
V p = primary voltage 

Ypri = admittance looking into the primary 

Y = admittance looking into the secondary 

sec 

Rp = primary load 

R s = L-L secondary load 

L * leakage reactance (referred to secondary) 

R = winding resistance (referred to secondary) 

w 

C - stray capacity (primarily bushing) 

Using this model, the equations for the quantities of interest 
become : 

v p N LR S + R w + jwLj 

Y . - jwC + —x 

P ri N Z [R c + R w + jwL] 

5 W 
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1 


JwC + 1 


= 

V s N [R w + jwL+-g-I 




rr^ 


N ( jwC + 4-) 


If jwC << y- and -E- «^jR^ + (jwL) 2 then 
k d it V w 


V p R p 

v; = W tR w + m3 


v sec * C R w + j" L + ^ —— 7 ^ 


N ( 1 /R p + jwC) 


If jwC << I— and 
R P 




+ (jwL/ then 


Y sec * (R w + J" L > 


•1 


The unknowns in these equations are R w , L and C. To check the 

validity of the model, it was necessary to determine values for these 

unknowns and thus was done as follows. 

To find C for the 25 kva GE transformer, we take the curve of 

Fig. 25 (~— ), which is the admittance looking into the primary with 
Z 11 

the secondary open-circuited, determine jwC for any frequency 

and then determine C at that frequency. If we believe the model, then 
this value of C will represent the C of the model and for this trans- 
former the value of C, measured at 5 kHz, is determined as follows: 
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/ 

= y. = -88 .5db X85 0 9 5 kHz 

z n in ‘ 

1. - 88 ,5db /85° 

i n 

Z. = 2319 - j 26506 
1 n 

C = 1.2 nf 

The small resistive component has not been identified, but may 
be attributed to measurement error. The value of C for the 25 kva 
Colt transformer can be similarly determined from Fig. 20 and is 
found to be 1.5 nf . 

We would expect the value of C measured at any other frequency 
from Fig. 25 to remain the same as that for 25 kHz and this is the 
case within reasonable limits. Figure 25 can be seen to be pretty 

much a curve of jwC versus frequency as shown in the jwC plot of 
Fig. 11. 

To determine the values for and L for the 25 kva GE trans- 
former, we take the curve of Fig. 17, which is the admittance looking 
into the secondary (L-L) with a 10 ohm primary resistance, and deter- 
mine Y at a given frequency. From this value of Y cor we can deter- 
sec ■ sec 

mine R, and L as follows: 
w 

Y * -21.7db/84° @ 25 kHz 
sec * 

Z co „ = 21.7db A 84° 

« 12.16 /84° * 1 .27 + j 12.10 
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and then 

R w - 1 .27n ) measured at 25 kHz 
L = 77 

Ideally, this measurement should be made with the primary short- 
circuited, however, the 10 ohm primary resistor becomes insignificant 
because 10 ohms « and, when transferred to the secondary, 10 ohms 
divided by N 2 (N 2 = 1000) is much less than jwL. 

Determining R and L at other frequencies from the curves of 

W 

Fig. 17 yield the following results: 


Frequency (kHz) 

R w (ohms) 

jwL (ohms) 

L (uh) 

5 

0.14 

2.72 

86 

10 

0.3 

5.29 

83 

25 

1.27 

12.10 

77 

50 

2.76 

23.8 

76 

100 

11.90 

58.0 

82 


Similar determinations can be made for the Colt 25 kva trans- 
former from the curves of Fig. 12, where Rp = 0, as shown below: 


Frequency (kHz) 

R w (ohms) 

jwL (ohms) 

l (uh) 

5 

0.23 

2.59 

82 

10 

0.68 

4.85 

77 

25 

2.21 

10.87 

69 

50 

4.75 

20.59 

66 

100 

21 .52 

51 .95 

83 
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With regard to the GE transformer, several observations can 

be made about both R and L. R is significantly less than jwL up 

w w 

to about 50 kHz and behaves, to a first approximation, according to 
the equation 



1.4 


where f = frequency in Hz 

R = R (ohms) measured at 5 kHz 
w w 


up to about 50 kHz. Computing R versus frequency using the above 

W 

equation yields the following result: 


Frequency (kHz) 

R (ohms) 
w 

5 

0.14 

10 

0.38 

25 

1.37 

50 

3.61 

100 

9.52 

It can be seen that the maximum apparent change in L with fre^ 

quency is less than 19% or less than 2 db, for either transformer. 

The Colt transformer behaves 

similarly to the GE except that 

its R. values are about twice that 
w 

for the GE which makes them more 

significant when compared to jwL. 

Computing R w versus frequency us 

the above equation yields the following results: 

Frequency (kHz) 

R (ohms) 

5 

0.23 

10 

0.62 

25 

2.25 

50 * 

5.93 

100 

15.64 
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A model of a distribution transformer has been developed and 
equations derived for the characteristics of interest as a function 
of frequency. The unknown quantities in the equations have been 
measured [for 2 transformers] and now it is appropriate to determine 
how well the model can predict measured performance. A computer 

plotting program was developed to accept R w> L and C as inputs and plot 

V V 

the equations previous given of , Y . , _j) and Y , as a function 

V P ' ' y sec 

of frequency, and to match the abscissa and ordinate dimensions as 
generated by the x-y plotter. 

Figures 28 through 31 show these curves for the GE 25 kva trans- 
former and are to be compared with the measured performance of Figs .14- 
17. The value used to generate the curves of Figs. 28-31 are as follows: 

R = 0.14ft @ 5000 Hz 

w 

L = 86 yh 
C = 1.2 uf 

The equation used here to determine R as a function of frequency 

W 

was implemented in the plot routine as 

f !- 4 
R w = 20,000 

instead of 

R ws f K3 
R w = .14 20,000 

For the 25 kva GE transformer, this makes a trivial difference 
since R is nearly ten times less than jwL at any frequency up to 50 kHz. 

W 
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A comparison of the computed curves with the measured curves 
shows the following: 

v s 

1. Both curves of — show agreement within 1 db up to 

P 

50 kHz . Incidently, the computed curve falls exactly 
on the 3 points computed from the algebraic model. 

2. For Y pr ^ , the resonance occurs at essentially the 
same frequency; however, the "Q" of the computed 
curve appears to be somewhat higher'and the maximum 
variation in magnitude, up to 50 kHz, is about 5 db. 

In this case, the algebraic model seems to more closely 

follow the measured curve than the computed curve. 

V 

3. The curves for track within 1 db up to 50 kHz. 

4. Y sec also shows good agreement between the computed and 
measured curves . 

A similar comparison can be made between the measured [Figs. 10- 
13] and computed curves [Figs. 32-35] for the 25 kva Colt transformer 
and the above statements, in general, appear to be applicable here. 

In conclusion, we believe that the simple "RLC" plus ideal 
transformer" model that we have developed, describes the required 
performance characteristics of distribution transformers , up to 50 kHz, 
with sufficient accuracy for the intent of this program. 
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Measured and Computed Performance as a Function of Load Changes 


Figure 36 shows the measured performance of rr- of the 25 kva 

V S 

Colt transformer when R p * lOOf. and R p * 1000ft. Figure 12 shows the 


measured performance when R p * lOf. and R p * 100ft. These can be compared 

V p 

with the computed performance of y— at R p * 1 Ofi [Fig. 34] and at 

s 

R p * 100n [Fig. 37]. The correspondence between measured and 'computed 


is reasonably good. 

Figure 13 "shows the measured performance of ^ sec * again for the 

25 kva Colt, as a function of R p and thus can be compared to the computed 

performance of Fig. 38 where R p * 100ft and Fig. 35 where R p s 10ft. 

Figure 32 shows, for the 25 kva Colt, the computed performance of 
V s 

y— when R<. = 4ft and Figs. 39 and 40 show this when R<- * 2ft and 8ft, 

V P b V S 

respectively . The measured performance of tj — when R,. s 4ft is shown 

V P b 

In Fig. 10. 

Figure 11 shows the measured Y . for the Colt transformer when 

pn 

R<- = 4ft and Fig. 33 shows the computed curve for this condition. 

Figures 41 and 42 show computed Y pr ^ for R^ = 2ft and 8ft, respectively. 

Again, the agreement in measured and computed performance up to 
50 kHz is reasonable. 
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Extended Frequency Measurement 

V 

Figure 43 shows the measured performance r the 25 kva 
GE transformer, from 5 kHz to 500 kHz for different values of R p . 
It can be seen that, up to 50 kHz, these curves agree closely with 
the corresponding measured curves of Fig. 16. It can also be seen 
that the performance above 100 kHz indicates significantly more 
complex behavior than described by our simple model. 
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Three Phase Transformer Measurements 

Some measurements have been made of 3 phase transformer banks 
on earlier programs and that Information is Included here only to show 
the kinds of things that have been done to date and perhaps to suggest 
further work that should be done in this area. 

In February of 1980, measurements were made on a number of 3 phase 
grounded wye to grounded wye pad mount transformers in the following 
manner. A resistor load was placed on the transformer secondary which 
consisted of 20% delta load and 80% wye load and with the total load 
being variable from no load to full load for each kva rating. Measure- 
ments were made at 4 frequencies from about 5-10 kHz and consisted of 
a Y.j i measurement and a i measurement which are defined according 
to the following diagram. 



where I A is the current measured going into 
the box and V A is the applied voltage measured 
to ground. 



rj— wi th A & C grounded 
V B 
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The results of these measurements are shown In Fig. 44 for a 
500 kva GE pad mount. 

Similar measurements were made on a 3 phase bank connected also 
grounded wye to grounded wye and consisting of three 75 kva dis- 
tribution transformers as shown in Fig. 45. Figure 46 shows the 
results of those measurements. 

The same kind of measurements were made on a 3 phase bank con- 
sisting of two 25 kva transformers in an open delta configuration as 
shown In Fig. 47. Figure 48 shows the results of those measurements. 

Figure 49 shows Y pri - for the 25 kva GE transformer previously 
mentioned when its secondary is connected In parallel for 3 phase wye- 
wye bank configuration. This is shown for various secondary loads. 
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Data Bank 


For purposes of feeder modeling, it is necessary to obtain the R , 

w 

L, and C values of all transformers connected to that feeder so that, 

in conjunction with the transformer model just aeveloped, the effect on 

the feeder of these transformers can be determined. Several ways of 

obtaining R , L, and C will be discussed, 
w 

The obvious and most accurate means would be to measure each 


1 

*17 


and Y sc as 


transformer. This would involve the measurement of 

previously described where is the admittance looking into the primary 

*11 

with the secondary open circuited. The value of C can be determined 
from y— . The measurement of Y $c consists of measuring the admittance 


'll 


looking into the secondary with the primary short circuited. The values 

of L and R can be determined from this measurement, 
w 

Another method is to take a statistical average of R ,, L. and C 

w 

for each KVA rating based on a large number of measurements of the type 
just described. 

A third method is to determine R and L from the nameplate data 

w 

and to obtain C from a statistical average. The nameplate of each 
transformer contains a Z value, for example Z = 1.9%, which represents 
the impedance at 60 Hz when looking into the primary, as a percentage 
of the full load rating for that transformer. It is not clear, in 
typical feeder modeling applications, that the nameplate data of in- 
service transformers is always available without actually looking at 
the nameplate which, in most cases, is rather difficult. 


As a part of this program, a number of transformers were measured 

reqard to R , L, and C using the y— and Y. . technique previously 
3 w 1 sc 

described. Table 1 shows a typical portion of that data. Table 2 is 


wi th 


a summation of the L measurements for different KVA ratings and 
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TABLE 2 

L(mh) Measured at 5010 Hz 
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Standard Deviation 
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manufacturers. Table 3 Is a summation of the C measurements, again 

as a function of KVA and manufacturer and Table 4 is the same for R . 

w 

In Table 3, the values under the "measured" column were obtained 
from a measurement of and the values under the "Bridge" column 

ni 

were obtained using an ECO Model 100 capacitance bridge from one primary 
terminal to the neutral connection. 

A comparison of C values obtair.td from both actual measurement 
and bridge measurement on the same 25 KVA transformers is shown below: 

1 


Manufacturer 

z n 

ECO Bridge 

Colt 

1 .47 nf 

2.1 nf 

McGraw-Edison 

1.25 nf 

1 .24 nf 

II II 

1.25 nf 

1 .25 nf 


The agreement between the 2 measurement techniques is good for 
the two McGraw-Edison' s , however, the bridge reading for the Colt is 
about 43% higher. 

Figures 50-52 show statistical values for L, R , and C as a 
function of KVA as obtained from Tables 2, 3 and 4. 

As previously mentioned, another method for obtaining R w and L 
is to use the Z percentage value stamped on the nameplate of all dis- 
tribution transformers. 

Table 5 shows R and X. (at 5010 Hz) obtained from nameplate data 

W L 

of all the transformers measured and compares the accuracy of this 
method with that for the statistical determination. 

Looking at Table 5, the Z column, in percent, is nameplate data 
and represents that part of the impedance, at 60 Hz, seen when looking 
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TABLE 3 

C (nf) Measured at 5010 Hz 



Standard Deviation 0.10 0.25 0.22 0.88 



















TABLE 4 

(ohms) Measured at 5010 Hz* 
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into the primary of the transformer, at rated load, which is contri- 
buted by the transformer itself. 

As an example, take the 1st transformer of Table 5, a 10 KVA 
Allis-Chalmers transformer with a nameplate Z value of 1.5%. To get 
this Z value in terms of 5000 Hz, which are the values listed in 
column 1, use the formula: 


Z NP-5000 


.015 ( 


240‘ 

10,000 


) ( 


5000 v 
60 ' 


Columns 3, 4, and 5 show the actual measured values of R , X 

w l 

and Z-for each transformer at 5000 Hz. A comparison of columns 1 and 
5 shows that the Z value obtained from the nameplate is larger than 
the measured value. No attempt is made here to explain this difference. 
Column 6 is a ratio of measured Z to the nameplate Z and averaging this 
ratio for all the 10 KVA transformers led to the .55 value shown in 
column 6. The nameplate Z's were then multiplied by this .55 factor 

I 

and column 2 was generated and called Z^p. 

Obviously Z^ p consists of both R w and but it is not obvious 

•n what proportion. Again, the measured ratio of X to R was computed 

w 

and this is shown in column 7. An average of the ratio, for 10 KVA 
transformers, was found to be 14.46 and this means that, at 5000 Hz, 
the R w is an insignificant part of Z and therefore the following is 
assumed to be true: 


X^ = Z^ p at 5000 Hz. 

Column 8 represents the difference between the value of X^ 
determined from the nameplate and the measured value for each transformer. 
Column 9 shows the difference between the value of X|_ obtained by 
averaging all the measured values of Xj_ and the measured value itself. 


A comparison of columns 8 and 9 shows that the average difference is 
less for the nameplate method but perhaps not significantly. 

While R w is not significant at 5000 Hz, it does become signifi- 
cant at higher frequencies and therefore, if nameplate data is to be 
used, a determination of R.. must be made at 5000 Hz so that it can be 
used in the previously mentioned formula 

R' -1.4 

D - f ) 

K w .14 '20,000- 

where R* - R (in ohms) measured at 5000 Hz 

w w 

To obtain R * , the Z obtained from the nameplate was divided by 
X L 

the averaged ratio of of column 7 and these values of R' are shown 

R w W 

in column 10. 

Column 11 shows the difference between the R, , values of column 10 

w 

and the measured R values, of column 4. Column 12 shows the di "erence 
w 

between the measured values of column 4 and the averaged values of 
column 4. A comparison of columns 11 and 12 show minimal difference 
between the 2 methods. 

It appears that taking a statistical average for the values of 
R and L does not appreciably degrade the results as compared to the 

W 

nameplate method and it is probably significantly easier to achieve. 

It is recommended that more statistical data be taken on trans- 
formers, especially with regard to L and C. 


TABLE 5 


A Comparison of Measured, Statistical and 
Nameplate Values for X L and at 5000 Hi 
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GE 
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Table 5 (cont'd) 

A Comparison of Measured, Statistical and 
Nameplate Values for and R w at 5000 Hz 
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Introduction 


Power distribution systems are increasingly being used for 
communications and control using power line carrier technology. For 
applications like meter reading, it is important to have a two-way 
communication path. 

Often, due to the length of sections of line involved, and due 
to branching, some points on a distribution feeder are difficult to 
access. This may occur due to standing wave nulls in the voltage profile. 
In such cases, it may be desirable to compensate sections of distribution 
feeders so as to improve reliability of communications to and from such 
points. In the present report we investigate the mechanisms of signal 
propagation on distribution systems and develop procedures useful in 
determining locations and types of compensations which can be used to 
enhance communication reliability on such systems. 

Inasmuch as the compensation approach we use is procedure oriented, 
it is useful to obtain an understanding of the factors which determine 
propagation transfer functions on networks having structures similar to 
those of distribution systems; distributed tree networks. To that end, 
we first consider propagation on a lossless section of two-wire line 
having a transformer on it and look at the factors ^hich determine pro- 
pagation to and from the transformer secondary, taking into consideration 
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typical properties of the distribution transformers. Following this, 
we develop compensation procedures on such a line. The results of this 
analysis are then extended to distributed tree networks of two-wire 
lines and subsequently to mul ti -conductor systems. In order to simplify 
analysis, we introduce appropriate approximations. These are motivated 
and justified from practical considerations and from experience gained 
in working with distribution system computer models and field experi- 
mental measurements. 

A. Simple, Lossless Two-wire Line with Transformer 

Figure 1 shows a lossless two-wire transmission line of length 
i , with a distribution transformer connected at a distance d from the 
input end. The transformer primary input impedance is assumed to be 
large in comparison with the characteristic impedance of the line 
Zo (* 1/Yo). Typical values of line characteristic impedances range 
from 100 ohms to 400 ohms, while primary side input impedances of 
distribution transformers range from 10 Kilo-ohms to 100 Kilo-ohms. 



In order to obtain the voltage transfer, one may proceed as 
follows. Let the secondary terminals of the distribution transformer 
be open circuited. Let a voltage vector of value Vi be impressed at 
the input. 

The voltage V 3 , developed at the secondary terminal of the trans- 
former, will be equal to the product of the voltage V 2 developed at the 
primary terminals of the distribution transformer and the voltage ratio 
of the transformer. The current , at the input end, will be the 
product of the input admittance Y] and the impressed voltage Vi . Inas- 
much as the individual transformer primary input impedances are large 
and do not significantly affect the propagation on the line, the voltage 
V£ at the transformer primary terminals will be determined principally, 
by the position of the transformer on the section. The voltage transfer 
across the transformer will be independent of the position of the trans- 
former. 

With the receiving end open circuited, if we define V R as the 
receiving end voltage, we can express 

V 2 = V R Cos 3 U-d) 

V 1 = V R Cos 3 A 

where d = distance from source to DT 

Similarly, if we express the impedances looking towards the receiving end, 
we can write 

Z 2 = jZo Cot 3 U-d) 

= jZo Cot 3 & 

Then, ^2 _ Cos 3 U-d ) 

V.j “ Cos 3 a 
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Now, Cos 2 x = - °~ X j 
Cosec^x 


Cot 2 x 7 
1 + Cot x 



Cos 2 B U-d) J Izf Z ° + Zl 

Cos 2 6H ^ Z 1 / Z 0 2 + Z 2 2 


Figure 1 shows a plot of Impedance magnitude normalized to Z 0 vs. 
voltage V 2 normalized to Vr # 



Figure 2 

The relationship is monotonic, indicating that high voltage 
transfer corresponds to a high impedance at load point and voltage 
nulls correspond to impedance nulls. The essence of any compensation 
scheme will be to eliminate voltage nulls. The above analysis indica- 
tes that, at least in the case of a single line section, this can be 
accomplished by the elimination of impedance nulls. 


\ 


t 


i 
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Reverse Path Considerations in a Two-Wire System 

In order to understand the factors which influence the reverse 
path propagation, we make use of reciprocity theorems which state that 
the trans-impedances and trans^imi ttances of a passive, bi -I inear 
networksare equal. For a two-port, specifically, the theorem states 
that Z 12 » Z 21 and Y^ 2 = Y 21 . 

The voltage ratio of the outbound path can be expressed as 



where I-jis the current flowing into the network at the input. The 
second term is the input admittance of the network. If we define Z 3 -j 
as the trans-impedance, then the above can be expressed as 



where Y^ is the input admittance at point V-| 

The reverse path voltage ratio can be expressed as 



where Z-j 3 is the trans-impedance from the secondary side and Y^ is the 
secondary circuit input admittance. Using the reciprocity relationship, 
we can express the reverse path voltage ratio as 

B R ■ e F • V Y 1 

Before proceeding further, we wish to remark that the above 
relationship is valid regardless of the complexity of the network 
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intravem'ng between the source point and the transformer secondary 

encountered. The above relation is also valid for lines with losses, 

and is not restricted to lossless lines, and can be extended to n and m terminal 
two ports. 

The above equation indicates that the reverse path voltage gain 
will be highest when Vp and Y 3 are high, arid Y-j is a minimum. The value 
of Y.j to be used is the value of the input admittance in parallel with 
the source admittance. 

Due to the presence of a relatively high leakage reactance at 
power line carrier frequencies, the secondary input admittance Y„ is 
relatively independent of the primary side admittance loading of the 
transformer. 

Of the two remaining factors, Y^ depends upon the admittance of 
the outbound source (or the inbound receiver admittance) and upon the 
line admittance at its location. Vp, as shown earlier, is related to 

the outbound line impedance at the point where the transformer is located. 

Thus, in order to obtain a high reverse path voltage gain from 
all points of the simple line section, it is necessary to have a high 
outbound line impedance at the different transformer location and to 
locate the receiver at a point of standing wave maxima and also have 
a high receiver impedance. 

Siqnal-to-Noise Ratio Considerations on a Simple Line 

Consider a simple line section with several transformers on it. 

The transformers are assumed to be randomly placed on the line section. 
Furthermore, it is assumed that noise originates on each transformer 
secondary with a mean value of v^ and a standard deviation of v a . 
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The noise amplitude distribution is assumed Gaussian. We wish to 
investigate the effects of line VSWR on worst case signal-to-noise 
ratio, of inbound signals. 

In the following, we will assume that the secondary input 
admittance Y^ is independent of transformer location. As before, we 
assume that the transformer does not significantly load the line. 

For a detector located at the point Y.| on the primary, the reverse path 
voltage transfer ratio, for the simple section, will be expressed by 


•k ’ S F V Y 1 

Bp, the forward path voltage transfer, will be a product of the 
transformer primary to secondary voltage ratio and of the voltage 
transfer ratio to the primary of the transformer. For a simple line 
segment, if K is the VSWR on the line, for any given primary receiver 
for inbound communications, the signal to noise ratio can be estimated 
as follows. 

Let N be the number of transformers, let the maximum reverse path 
gain be G mav , and let |pj be the magnitude of the reflection coefficient. 

fliaX 

The range of reverse path gains, G, will vary between 


(1- 

p_ 

2 

: 1 + 

p 



< G < G 
— — max 


Since |p| is related to the VSWR K by 


K-l 

K+l 


we can express the above relation as 



< G < 


G 


max* 
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If V g is the amplitude of the secondary signal (assume constant 

on every secondary), the worst case signal amplitude in the receiver 

V • _ will be G. • V 
mi n mi n s 


i .e. , 


V G 

_ s max 


min 


If Y Qr is the real part of the receiver admittance, the minimum signal 

2 

P° wer P smin wi11 be V m in V 1 ' 6 '* 

V 2 G 2 

p - s max Y 
smin k 2 or 

In order to obtain the mean noise level, we consider the 
relationships represented by the phaser diagram shown below. 



Figure 3 

OA 

PA represents the maximum gain G max , the ratio gp is the 
magnitude of the reflection coefficient Jpj. PC represents the minimum 
gain. PB is the gain to a general point. 
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(PB) 2 = (OP) 2 + (OB) 2 - 2 (OP) (OB) cos e. 


2 

The average value of (PB) will be 


2u 

27 \ (PB ) 2 de 


2tt 


2ir 


(OP) 2 de + ^ 


2tt 

( 


(OB) 2 de 


J 

0 


1_ 

IT 


2tt 


(OP) (OB) Cos 8 de 


i .e. , 


<(pb)^> =«op) 2 >+4ob)^= ^ 


pI 2 (PA ) 2 


(1 + |o| ? " (1 + 1 p i 2 


* < pa > 2 

(t + Ipjf 

where < > denote the expected value. 

The mean noise power will be: 

p = 1 kL g 2 n( v ^+v 2 )Y 

N (l + |pjP maX m ^ ° r 

■ ^ G .a X 2 » + V 2 > Y or 


Thus, the ratio of minimum signal power to mean noise power 


will be 


s min „ 
P N ' 


2 

Tk 2 " 


tt 


m 


+ v 
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The ratio 


2 “2 

v m v cr 


is the ratio of the signal to mean 


noise powers on each secondary. 
Thus 


J min 


receiver 


N (K 2 + 1) 




secondary 


We should make the following observations. The above result 
indicates that the worst case signal-to-noise ratio will depend upon 
the VSWR and will be poor for high VSWR. The dependence on VSWR is, 
in reality, only a dependence on the ratio of maximum to minimum gains. 
For an uncompensated line, this ratio is the voltage standing wave 
ratio. Thus, in order to operate the line so as to maximize the worst 
case signal to noise ratio, it *■% .-^cessary to reduce the ratio of 
maximum to minimum gain on the line. 

The analysis suggests that the signal to noise ratio, as far as 
noise generated from the secondaries, is independent of the location of 
the receiver and its impedance level. On the other hand, location of 
the receiver and its impedance level will affect the received level of 
noise generated on the primary. 

The analysis of this section is not limited to simple lines, and 
can be extended to branched lines. In the latter case, the concepts of 
minimum gain of signal to average gain of secondary noise will yield 
worst case signal to noise ratios. 
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Compensation of Simple Lines 


We first investigate the compensation of simple lines so as to 
provide guidelines for approaches to be used in branched situations. 

The purpose of compensation is to provide maximum reverse path signal 
to noise ratios. Ideally, the elements used to compensate the line 
should draw minimal amounts of power from the signal, and from the 
power line. 

As noted earlier, the worst case signal to noise ratio depends 
upon the ratio of maximum to minimum gains and for maximizing S/N 
ratios it is desirable to reduce this ratio. In a simple line, this can 
be achieved, in theory at least, by terminating the section with its 
characteristic impedance. There are two problems with this approach. 

On one hand, the characteristic impedance of most line sections is 
resistive (dissipative) and of a relatively low value, ranging from 
15fi (for underground cables) to approximately 400fl for overhead lines. 

In order to overcome primary generated noise, this entails the use of 
large secondary signal transmitter power. On the other hand, impedance 
mismatches between overhead and underground sections will create standing 
waves with large VSWR. 

Experimental evidence obtained during field verification of the 
DIFNAP model indicate that as far as signal propagation on the primary 
is concerned, the distribution line behaves like a lightly loaded trans- 
mission line. The line parameters we have computed indicate a propa- 
gation constant tangent angle (tan-1 0 f ) of less than 2° for overhead lines in 
creasing to about 6° for underground cables. Even with loading, the 
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effective propagation constant tangent angle usually is less than 10°. 
A line terminated by its characteristic impedance will require from 
10 to 20 times the signal power to establish the same voltage level 
as an unterminated line. Some standing waves will also be created at 
branch points. 

Another approach is to employ reactive shunt compensations. In 
order to understand this method of compensation, let us consider a 
short length of 2 wire lossless line open circuited at its end. The 
length & is less than a quarter wave-length. 

t = — t - ==d 

r-> 


Figure 4 

The input admittance will be 
Y, n * jY 0 tan 8* 

where Y' is the characteristic admittance of the line section and 3 is 
the propagation constant, at the frequency of operation. Let a lossless 
compensation of admittance jY be placed at the midpoint of the line 
section. The admittance of the line to the right of the compensation 
will be 

Y r - jY o tan 

The sum of this admittance, together with that of the compensation, 
will serve as the load on the left hand section of length lit ^ and the 
overall input admittance can be calculated using transmission line 
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formulae to be 


Y 2 cos + j Y q sin |i 

Y*T (with compensation) = Y 

ln 0 Y o cos |* + j Y, sin P 


Substituting for Y^ = j[Y c +Y q tan ] 
we obtain 

. jY cos P + j 2 Y rt sin P 

w C _ y C 2 Q Z_ 

in o w ... 32. » • B i 


Y cos ^ - Y sin ^ - Y rt sin ^ tan P 
o 2 c 2 o 2 2 


B2 


Dividing numerator and denumerator by cos f- o 


we have 
,c 


32. 


YlT = jY ' Y c + 2Y o * an 2 
in o TbF” 


Y o [l-tan‘ - Y c tan ££ 


If Y , —2 Y_ tan P then 

C O 2 


Y. c 

in 


=JV o - 2 Y 0 tan + 2 Y 0 tan ^ 
V 0 C 1 + tan 2 S|] 


= 0 


The magnitude of the maximum admittance along the line section will be 

Y t an p # What has been achieved is that the admittance on the line 
o 2 

looking toward the receiving end will vary from zero to a (+) maximum 

0 

magnitude of Y c just at j + t and back to zero at the sending end. 

By dividing a given length of line L into n equal line sections, 
each shorter than and placing a compensation at the center of each 
section, it is possible to make the maximum admittance to be 

Y rt tan < 0.41 Y . 
o 2n o 
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The ratio 


max 


mm 


4 which is equal to the ratio of maximum to 
SL 


minimum gains will be Sec ^ < 1.083. 

If each section was approximately a quarter wave long, the ratio 
of maximum to minimum gains would be 1.414, with this approach. The 
plot of Fig. 4A shows tlie relation between the electrical length section 

and the ratio of maximum to minimum qains. 

It should be remarked that at one frequency, the ratio K of the 
maximum to minimum reverse path gains will hold regardless of the 
overall length of line as long as the individual compensated sections 
have the corresponding maximum electrical length. As shown below, this 
procedure can be generalized to tree configured networks. 


Generalization of Compensation Procedure to lossless Tree Networks 

Consider a lossless 2 wire tree network as shown in Fig. 5. 
Without loss of generality, we assume that the receiver is located in 
the root section of the tree. Transformers are distributed in a random 
fashion throughout the network. As before, we assume the generation 
of noise with a Gaussian Amplitude distribution with an R.M.S. mean 

i, L 

value v and an R.M.S. standard deviation v . Let T. be the i number 
m a i 

of the total N transformers on the system. 
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Let Gi denote the reverse path gain from the i^ transformer to 
the receiver. As noted earlier, this gain will be directly related to 
the forward path gain, between R and Ti, along the primary circuit. If 
J1 , J2, — JP denote the junction points on the path from the receiver 
to Ti, then the forward path voltage gain Vp^ from the receiver to Ti 
can be expressed as 


V 


Fi 







In general, because of branching, many terms in the above product 

will have a magnitude less than 1. A few of the terms, especially 

on sections at the end of the tree, may have magnitudes greater than 1, 

The worst case gain, Vp will result from Vy. being at a voltage 

mi n 

minimum of a VSWR and a low combination of transfers to the branch point 
at the source end of the branch on which the corresponding transformer 
is located. The worst case reverse path gain G . will correspond to 
such a transformer. 

The worst case signal to noise power will be 

„2 


smm 

> 

noise 


G 2 . 

min 

N n 2 

z: G i 

i=1 


(V 2 + V 2 ) 
m a 


In order to obtain an optimum worst case signal to noise ratio, 
we have to maximize G - n in magnitude. This implies maximizing the 
forward primary level path gain. This can be accomplished by minimizing 
the maximum to minimum gain ratio on each branch, treated as an isolated 
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branch and also maximizing the voltage transfer across each branch. 

The former can be accomplished using the procedure outlined in the 
previous section, while the latter can be achieved by minimizing the 
magnitude of the input admittance to each branch. 

Extension to Two-Wire, Loaded, Lossy Lines 

In order to extend the above analysis to real life situations, it 

is important to realize that power distribution feeders, at the carrier 

frequencies, are low-loss transmission lines and the loading on them 

corresponds to a lightly loaded situation. As a result, they exhibit moderately 

high VSWR. The analytic treatment of the problem is relatively complex, 

even for a two-wire line using low loss approximations. 

Before we begin to extend the analysis to the lossy situation, it 

is instructive to observe the effect of compensation of lossless lines, 

analyzed earlier, on a Smith Chart. Figure 6 on the Smith Chart shows 

the impedance paths along the uncompensated and compensated paths for 

a line of electrical length 0.18 wavelengths. The maximum to minimum 

voltage ratios for the uncompensated and compensated cases can be 

determined by using the cosines of the appropriate electrical angles. 

For the case illustrated, the angles for the uncompensated case are 0° 

at the beginning of the section and 65° at the end of the section. 

The V„, to V ratio for the uncompensated would be 1: 0.422. For 

the corresponding compensated case, the angles would be 0° and + 32.5° 

and the to \L„- ratio would be 1 :0.843. A 6 db improvement in 
max min 

worst case signal to noise ratios will result, according to the 
previous analysis. 
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In extending the analysis to the lossy, loaded situation, we make 
the following observations: 

a) The path along the lossy line will represent a spiraling 
circular path along the Smith Chart. 

b) The lumped transformer load admittances will create small 
discontinuities in the path. The aggregate effect of these 
discontinuities will be to change the nature of the spiral in 
a small way. 

c) Because the spiraling (which occurs as a result of the lossy 
medium), the input admittance of the compensated section 
cannot be made zero by the use of a lossless compensation. 

In order to make the input admittance zero, the compensation 
would have to add the energy consumed by the loss of the 
network. 

d) The implication of the previous observation is that in 
subsequent sections the section load admittance will be non- 
zero. At a single frequency, it may be assumed that this load 
admittance will be real and will have a magnitude much smaller 
than that of the characteristic admittance of the line. 

In order to demonstrate the effect of loading with the 
Smith Chart, consider the loaded section shown in Figure 7. 


1 1 J 
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rr — i 
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Fig. 1.20 b 

Figure 8 j 

Profile of Admittance Along a Loaded, Lossy Section i 
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In Figure 7, represents the terminal load admittance which 
we assume is non-zero and almost resistive with a magnitude sub- 
stantially less than that of Y . The admittances of the loads connected 

(representing distribution transformers) have magnitudes much smaller 

« 

than that of Y . The plot on the Smith Chart of Figure 8 represents 
the impedance profile along the line section. The discontinuities repre- 
sent the effect of lumped admittances of the load. The diagram, since 
it is used to demonstrate the effects, has been deliberately exaggerated. 

Two observations need to be made regarding this diagram. The dis- 
continuities, since they represent passive elements, will always shift 
the curve towards the right, i.e., in the direction of increasing 
Re(Y). In the top half of the chart they represent a decrease in 
the magnitude of the reflection coefficient, i.e., they tend to push 
the spiral inwards. From the considerations of the previous section, 
it is clear that we will prefer to operate in this region of the Chart. 

The second observation which should be made is that the process of 
computing voltage transfer ratio in the case of the loaded line is 
somewhat complicated. As can be seen from the previous sections, 
control of voltage transfer ratios will be the principal objective Gf any 
compensation procedure. In order to explain this computation, we will 
demonstrate it using one discontinuity and exaggerate the effect of the 
discontinuity on the Smith Chart of Figure 9. The corresponding circuit 
is shown In Fig. 10. 



Figure 10 
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In the circuit of Fig. 10, represents the terminal admittance 

and Y. represents the loading admittance. The voltage ratio ^ is 

' V* v 2 

the same as the ratio of the vectors _A, Similarly, the voltage 

ratio ^2 is the same as the ratio lH. ' Thus, the voltage ratio li 
V 3 Vq v V 3 

will be J the product of the voltage ratios V A . . In the case 

VT V D 

of multiple load discontinuities, the products of the voltage 
transfer ratios across each section must be taken. 

This procedure can become complex, even on a Smith Chart, and the 
problem can quickly become intractable. In order to obtain a better 
handle on the situation, we use the following approach. Consider a 
section of lossy loaded line, with a compensation to be placed somewhere 
near the middle of the line. The two sections of line, on either side of 
the compensation, can be looked upon as lumped sections having a set of 
circuit parameters, as shown below. 



For a given load Y^, one can always determine the input admit- 

D 

tance Y in of the right hand section. With a lightly loaded short 
line and a terminal load Y L which is real and has an admittance 
magnitude much less than the characteristic admittance of the 
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Figure 9 

Demonstration of Computation of Voltage Transfer Ratio 
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D 

transmission line, Y in will be largely capacitive admittance with a 
small real part. 

D 

We can express Y- n as G-j + jS^ where G^ and are both real and 

G^<<Sp The load on the left hand section can be expressed as G^+j S-|-jY c 

where Y is the admittance of the inductive, non-dissipative, compensating 
0 

element. 

The input admittance Y^ n to the entire section can be expressed 
as 




iQ 2 

(y 22 £ + G-, + jS 1 -jY c ) 


where y. .1 are the short-circuit admittance parameters of the left 
^ J 

hand section. It should be clear that in the above expression, Y Q , 

S-| and G-j are real while in general, the short-circuit admittance 
parameters are complex. It should also be clear that with a suitable 
choice of Y c » the imaginary part of Y^ n can be made zero. From the 
discussions of the previous sections and comparison with the lossless 
case, it should be clear that with a choice of Y c that makes Im Y^ n 
equal to zero, the magnitude of the admittance will be a minimal. 

Generalization of Compensation Procedure to Lossy Tree Networks 
The compensation procedure of lossy, loaded sections can be 
extended to lossy and lightly loaded tree sections. The extension is 
analogous to the extension of the compensation procedure for lossless 
section to lossless tree networks. The compensation of each major 
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branch at the tree permits control of the ratio of maximum to minimum 
voltages on the branch. The input admittance to each major branch is 
minimized in magnitude and is made nearly real at and near the opera- 
ting frequency. In addition, the voltage transfer ratio across each 
major branch is controlled to have a magnitude close to unity. As a 
consequence, the ratio of maximum to minimum voltages on the entire 
network is controlled by the procedure to approximately equal the 
ratio of maximum to minimum voltage on any one branch, and a corres- 
ponding control of signal to noise ratios is obtained. 

Compensation Over a Band of Frequencies 

The preceeding discussion has been expressly directed to single 
frequency or narrow-band operation. The passive compensation schemes 
proposed will provide control only for a narrow range of frequencies. 

The use of passive compensation can become ineffective at certain 
frequencies, as can be seen from the following example. 

Consider a two-wire lossless line operating over a frequency 
range f 1 to 2fj . .Let us assume that we wish to compensate this line 
and choose a procedure so that the compensators are placed so that the 
first one is placed at a distance of ^ from the open end at the highest 
frequency. Furthermore, we choose to make the admittance at midband 
frequency minimal at every quarter/fcftv_e section at the highest frequency 
from the end. The admittance values normalized to Y Q at various points 
on the compensated section are tabulated below at each of the frequencies. 
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Distance from open 
end in Wavelength 
at highest freq. 

Admittance normalized to Y 0 
f ] 1.5^ 2f ] 

Condition 

0 

0 

0 

0 


A /8 

j 0.4142 

j 0.6681 

jl 

before compensation 

ty8 

- j 1 .595 

-j 0.6681 

-j 0.004 

after compensation 

X/4 

-j 0.71 

j o 

+j 1 .05 


3V 8 

-j 0.218 

j 0.6681 

j 50 

before compensation 

3V8 

-j 2.226 

-j 0.6681 

j 49 

after compensation 


Clearly, at the higher frequency, an impedance null will be 
observed close to the second compensation. 

Thus, using simple inductive compensations at the mid-band 
frequency is not an adequate compensating mechanism to control VSWR for 
broadband operation. Some improvement in the above situation may be 
realized by the use of compensators whose admittance vs. frequency 
characteristics are different from those of simple inductors. 

Alternatively, we may ask: what type of compensator can be used 
so that the input admittance to a section, having a length £ is nearly 
zero at all frequencies in the band? In the lossless case, assuming 
that the compensator is placed in the middle of the section, the com- 
pensator should have an admittance vs. frequency characteristic of 

Y = -j2Y Q tan ^ over the band. If £ is small such that tan can 

be approximated by over the band of frequencies involved, an element 

having the characteristics of a negative capacitor may be used. Such 

elements can be generated using active circuits and may be used if 

proper coupling practices are followed and appropriate drain inductors 

are used to ensure system stability at the power frequency and at 
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harmonic frequencies. Sections having lengths approximately equal to 
j at the highest frequency may be compensated in this manner. 

Extension to Polyphase Lines 

The compensation procedures outlined above can be extended to 
polyphase lines. The compensator used will be six-element devices and 
computations will be based on complex matrix arithmetic. The procedures, 
which need further investigation are clearly beyond the scope of the 
present report. They will involve the definition of appropriate opti- 
mation criteria and the development of computer algorithms to deter- 
mine optimum selection and siting of compensators. In addition, 
sensitivity to load lisle parameter variation and network reconfigura- 
tions will have to be investigated. 

Conclusions 

In this report, we have investigated the problem of compensating 
distribution lines, with transformer loading, so as to improve their 
performance and communication media in two-way communication systems. 
After investigating the principal phenomena which are responsible for 
the degradation of the medium as a communication channel we have 
developed criteria which can be used to determine the performance of 
the channel. Using these criteria as a basis, we have evolved methods 
of compensating the Tines. The methods, first developed for compensa- 
ting lossless, simple lines, are extended to the compensation of loss- 
less tree networks. The procedures are then extended to lines with low 
losses and loads and to lossly, loaded tree networks operating at a 
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single frequency. Finally, the compensation procedure is extended 
to broadband operation and to polyphase networks. It is shown that 
by appropriate compensation schemes, lines can be conditioned to serve 
as reliable two-way communi cation media. 
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